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ABSTRACT 

In this paper we analyze previously published spectra with high signal-to-noise ratios of E and SO 
galaxies in the rich cluster CL1358-I-62 at z = 0.33, and introduce techniques for fitting stellar pop- 
ulation models to the data. These data and methods will be used further in a larger study of the 
evolution of absorption line strengths in intermediate redshift clusters. Here we focus on the 19 ellip- 
tical and lenticular galaxies with an homogeneous set of eight blue Lick/IDS indices. These early- type 
galaxies follow very narrow line strength-line width relations using Balmer and metal lines, indicating 
a high degree of uniformity in their formation and enrichment hi stories. We explore these histories us- 
ing recently published, six-parameter stellar population models (jThomas. Maraston. fc Bendeill2003at 
IXhomas. Maraston. & Korn 2004), and describe a novel approach for fitting these models differen- 
tially^ such that the largest sources of systematic error are avoided. The results of the model fitting are 
accurate relative measures of the stellar population parameters, with typical formal errors of < 0.1 dex. 
The best-fit models yield a mean of 1.2 per degree of freedom, indicating that the models provide 
good descriptions of the underlying stellar populations. We find: (1) no significant differences between 
the best-fit stellar population parameters of Es and SOs at fixed velocity dispersion; (2) the stellar 
populations of the Es and SOs are uniformly old, consistent with results previously published using the 
fundamental plane; (3) a significant correlation of [Z/H] with galaxy velocity dispersion, in a manner 
consistent with the observed B — V colors of the galaxies, and indicating that dust is not a significant 
contributor to the colors of early- type galaxies; (4) a possible, modest anti-correlation of [a/Fe] with 
velocity dispersion, with < 2a significance, and discrepant with the correlation inferred from data on 
nearby galaxies at the < 3cr level; and (5) a significant anti-correlation of [ck/N] with galaxy velocity 
dispersion, which we interpret as a correlation of nitrogen enhancement with mean metallicity. Nei- 
ther [a/C], nor [a/Ca] shows significant variation. While the differences between our conclusions and 
the current view of stellar populations may point to serious deficiencies, our deduced correlation of 
mean metallicity with velocity dispersion does reproduce the observed colors of the galaxies, as well 
as the slope of the local Mg-cr relation. Our tests indicate that the inferred population trends do 
describe real galaxies quite well, and matching our results with published data on nearby galaxies, we 
infer that the discrepancy stems largely from the historical treatment of broadening corrections to the 
narrow indices. The data also strongly indicate that secondary nitrogen is an important component 
in the chemistry of elliptical and lenticular galaxies. Taken together, these results reduce early-type 
galaxies in clusters to a family with one-parameter, velocity dispersion, greatly simplifying scenarios 
for their formation and evolution. More specifically, our data conclusively show that cluster SOs did 
not form their stars at significantly later epochs than cluster ellipticals of the same mass, and the 
presence of secondary nitrogen indicates that both Es and SOs formed from self-enriching progenitors, 
presumably with extended star- formation histories. 

Subject headings: galaxies: clusters: individual (CL1358-I-62), galaxies: stellar content, galaxies: 
elliptical and lenticular, galaxies: evolution 



1. INTRODUCTION 

Ever since the disco very that early-type gal axies fol- 
low scaling relations ijMorgan fc Mavall Il957j) . it was 
clear they had profound implications for cosmology as 
well as for galactic structure, formation, and evolution 

^ Based on observations obtained at the W. M. Keck Obser- 
vatory, which is operated jointly by the California Institute of 
Technology and the University of California. 
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(e.g. iMinkowskil Il962t iTerlevich et"all 119811 and sub- 
sequent literature). Perhaps the most well-studied of 
these scaling relations is the color-magnitude relation 
()Bauml[l959y) . and its interpretation as a c orrelation o f 
metal abundance with galaxy luminosity ijRoodl IlijBflD 
has survived to the present-day. Over the past several 
decades, numerous other scaling relations have been dis- 
covered and interpreted, such as L-a (Fabcr & Jacksor 
197^, rp-Ip-a (the funda mental plane; Drcssler_£^^n 
1987. iDiorgovski fc David Il987jl . Mg-cr fTerlevich et al. 
Il981|) . and other line str ength-line width relations, such 
as 'U S A- or H'^a-o' (e.g. iKuntschneil l200nt iKelson et alJ 
MM)- 

Several important breakthroughs have helped to con- 
strain the nature of early-type galaxy scaling relations. 
The clearest results have come from explicit tests for 



2 



Kelson et al. 



changes in the ages of stellar population with redshift 
along the sequence of early- types. By measuring the 
slopes of the color-magnitude relations in clusters to red- 
shifts of unity 'Stanford. Eiscnhardt. & Dickinson (1998) 
concluded that the relation originates largely from sys- 
tem atic variations in met allicity with galaxy mass (also 
see iBlakeslee et aT] l2006j) . Likewise, the fundamental 
plane of early-type galaxies has a slope that does not 
appear to evolve s ignificantly with redshift to z = 0.33 
(jKelson et al.l2 000c'). further evidence that cluster E/SOs 
are a family well-described by uniform ages and a mass- 
metallicity relation. At redshifts of z = 0.8 — 0.9 there 
are hints that th e slope of the fundarn ental plane may 
evolve modestly (jJoreensen et al.ir2006|) . though this re- 
sult appears to be inconsistent with the slope of the 
color-magnitude relation as measured by Blakeslee et al.l 
((2006|) and it remains to be seen how these data will be 
reconciled. 

There have been many efforts to observe 
the high-rcdshift line stren gt h-line width rela- 
tions (Zieglcr fc Bended Il997t Kelson et all l2001t 
Jorgensen et all 12005: Ba rr et alJ 120051 iMoran et all 



2005ir However, significant constraints have not been 
forthcoming because of the difficulty of obtaining 
sufficient signal-to-noise ratios. 

The situation is much better at low redshifts. Detailed 
analysis of line stren gths in low-redsh ift E/SOs have been 
very revealing since iGonzalej lll993D f and other s) broke 
the age-metallicity degeneracv l|Worthevlll99^ . Using 
line strengths many authors have shown that massive 
galax ies in clusters are uniformly old (e.g. J0rgcn senl 
I1997L IT999I: iTrager et al.il2n00al: lKuntsclinerll200Ci) . but 
some ambiguities do remain. Fo r example it ha,s long 
been unde rstood that the Mg6-g ijTerlevich et al.llT98lj? 
and Fe-cr l|Gonzaled Il993j) relations are not consistent 
with the hypothesis that they both arise solely from 
a correlation between the mean metallicity of the stel- 
lar populations and galaxy mass. The discrepancy has 
been interpreted using stellar populations models in 
which "a-enhancement, or [a/Fe], is allowed to vary sys- 
temat ically wi th velocity dispersion (see, e.g. j0rgcnsCT . 
Rp9; Tragcr et al. 2000a b; Wo rthev fc Collobertfl20da 
iThomas et al...2005y for details and ma ny important ref- 
erences . While iTrageT et al.l (^2000a^ reminded read- 
ers that this is a misnomer, because the a elements 
are not enhanced — Fe is simply under- abundant. We 
will, however, continue to use the term "a-enhancement" 
to be consistent with past literature on the sub- 
ject. After several decades of analysis, improved 
techniques of observation, and despite the immense 
amou nt of progress i n stellar p opulation models (Fabcr 
1972; Worthev 1994; Tragcr ct al. '200'0a: Schiavon et all 
2002: Thomas. Maraston. fc Ben der 2003a) , the mod- 



eling of absorpti on lines r e mains ambiguous and con- 
tradictory (e.g.. iWqrl^h ev' '1994f; "Trager et al.' '2000, ' 
[ Worthev fc CollobertI 12003: Schiavon. CaldwcU, & Rosj 



120041 IThomas. Maraston. & Bender 2003b, and others). 
Some h ave even questioned the "a-enhancements" alto- 
gether l)Proctor et alJl2004a|) . 

Over the long-term we hope that the physics of stellar 
atmospheres will be understood with sufficient detail to 
allow for the creation of synthetic stellar spectra for stars, 
over the full range of stellar masses and phases of stellar 
evolution, that comprise the observed spectra of galaxies. 



Presently, however, numerous uncertainties in the forma- 
tion of molecular lines, and incomplete line lists prevent 
one from generating completely synthetic spectra for stel- 
lar populations. As a result the features in the spectra 
of stars and galaxies cannot be fully modeled. Despite 
major progress in generating high resolution sp ectral en- 
ergy distributions of simple ste llar populations ijVazdekid 
Il999t iBruzual fc Charloti2003f) , we cannot fully solve for 
the age(s) and elemental abundances of populations in a 
galaxy by directly fitting synthetic spectra. 

Without the ability to directly model the spectra of 
galaxies, one must measure and model spectral indices. 
On e widely used system is the Lick/IDS system, created 
by iBurstein et aL l)1984|) (and subsequently revised by 
iTrager eT^LflTgg? ). To this day, the measurement and 
modeling of these line strengths remain the most use- 
ful means of assessing the bulk properties of the stellar 
populations in passively evolving galaxies. 

Over the last several years several groups 
have expanded our ability to model these in- 
dices with more parameters than age, metallicity, 
and g - enhancement. The models o f ITrager et alJ 
(l2000aD IThomas. Maraston. & BenTIeil l|2003a() . aiid 
Schiavon (200^ allow one to probe specific elemental 
abundances, such as nitrogen, carbon, and calcium. 
Such models only make predictions for passively evolv- 
ing stellar populations. Fortunately models of galaxy 
evolution involving star-formation are not needed in our 
study of massive cluster galaxies through intermediate 
redshifts, because such galaxies have been shown to 
be passively evolv ing (e.g., ,^lson et al. 1997, 2000(J 
iWuvts et all l2004|) . At higher redshifts, on-gomg 
star-formation may become increasingly important (e.g. 
lJuneau et alJl2005() , and more sophisticated models may 
then be required. 

In our survey of galaxies in rich clusters at intermediate 
redshifts, the sample of galaxies with the highest signal- 
to-noise spectroscopy over the widest range of galaxy lu- 
minosities is that from our fundamental plane survey of 
CL1358-I-62 at z = 0.33. Because of the depth and qual- 
ity of that sample, we adopt it as the reference sample to 
which the oth er clusters in our su rvey will subsequently 
be compared (^Kelson et al.ll2006bj) . 

Here we study the absorption lines of the 
E/SO g alaxies in that sample, us in g the mod- 
els of IThomas, Maraston, fc Bended l)2003all and 
iTh omas. MarastoiT'X'lCornl lj2004jl to explore not only 
ages and metallicities, but additional stellar popula- 
tion parameters, namely the relative abundances of 
nitrogen, carbon, and calcium, along with the mean 
a-enhancement. Because we are keenly interested 
in potential variations in the abundance ratios, and 
because the state-of-the art high-resolution models 
(jBruzual fc Chariot I2003|) do not reproduce, in detail, 
the strengths of many features in the spectra o f real 
galax ies (e.g. CN, Ca4227, Mg, H5, H7; Gallazzi eTall 
I2005j) . we prefer to derive stellar population parame- 
ters from an analysis of absorption line indices. The 
IBruzual fc CharlotI l|2003() SEDs are employed for other 
purposes and these are discussed below in the context of 
deriving broadening corrections to our data. 

Many of the details regarding the processing of the 
data, and subsequent corrections, are discuss ed else- 
where (iKeioTerS] llOOOB jKeionir^ How- 
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Fig. 1. — Color-magnitude diagram of the cluster in R and Veoe — ^814) taken from iKelson et all I2000bl) . All 1 94 confirmed cluster 
members in the HST imaging are shown. The red, blue, and green symbols represent those galaxies observed by IKel son et al. (2000b) for 
use in their study of the fundamental plane. Ellipticals are shown using red, filled circles. The SO galaxies are shown using blue open 
circles. Spiral galaxies are shown in green. The light gray points were not observed in the high-resolution study of Kelson ct al. 12000b). 
The red, blue, and green points are shown using two sizes. The larger ones mark those galaxies in the homogeneous sample (for which all 
eight of the indices used in the modeling are available; see (0. The E and SO (red and blue) gal axies in the figure are discussed in this 
paper while the full fundamental plane sample of CL1358-I-62 will be discussed in a subsequent paper. 



ever several of the key points are discussed below 
in ^ because of their crucial role in the analy- 
sis. We then describe our comparison of the ab- 
sorption h ne strengths of the E/SO ga l axies t o the 
models of 'Thomas. Maraston. fc B ended lj2003all and 
[Thomas. Maraston, & Korn (2004), in which we perform 
non-linear least-sq uares fits to each g alaxy's set of line 
strengths (see, also lProctor et al.l2004j) . but only for rela- 
tive differences in the stellar population parameters . This 
allows us to explicitly avoid potentially large systematic 
uncertainties in the direct comparison of the data and 
models. This methodology is used in the rest of the sur- 
vey and so is described in some detail. The resulting 
relative ages and patterns of fitted chemical abundances 
are described in ^Ql^nd then discussed further in fJS] Our 
conclusions are summarized in ^ While our findings 
do not depend on the cosmology, we use the cosmolog- 
ical parameters HO = 72 km/s/Mpc, ^Im = 0.27, and 
fl\ — 0.73, when such parameters are required (e.g. for 
corrections to the line strengths for any variation in the 
metric sizes of the apertures from which the spectra were 
obtained). 

2. THE DATA 

As stated earlier, the galaxies in CL1358-I-62 were orig- 
inally targeted as part of a survey to study the fun- 
damental plane in intermed i ate redshift clusters (pub- 
hshed by [ Kelson e t al1ll997t Ivan Dokku m et"anil998bt 
IKelson et al. 2000c: IWuvts et al. 2004). The selection of 
the CL1358+62 spectrosco pic sample, and data process- 
ing, were fully described in iKelson et alJ l|2000b|) . and we 
summarize it here. 



The field of CL1358+62 was targeted for extensive 
imaging with the HST WFPC2. A two-color mosaic of 
the cluster CL1358-I-62 covering 64n ' was constructed 
using 12 pointings. Of the s pectroscopically confirmed 
members (jPisher et al.l 11998V 194 fall within the field 
of view of the WFPC2 imaging. From this large cat- 
alog, we selected a sample for detailed stud y with the 
Low Resolution Imaging Spectrograph (LRIS: l0ke et alJ 
Il995|) at the W.M. Keck Observatory. The high resolu- 
tion of the HST imaging allowed us to derive structural 
parameters with the accuracy needed for the fundamen- 
tal plane ifKelson et alJl2000ald|) . 

The sample used in the fundamental plane analysis 
was randomly selected from a catalog of known members 
within the field of the HST mosaic down to i? < 21 mag. 
Because the selection did not rely on morphological infor- 
mation, the galaxies span a range of morphologies , from 
ellipticals to spirals. In Figure^ taken from .Kelson et alJ 
(2000b|), we show the color-magnitude diagram for the 
spectroscopically confirmed cluster members within the 
HST mosaic. At the redshift of the cluster, L* corre- 
sponds to i? ^ 20.1 mag (adopting t he M/Ly evolution 
as measured bv IKelson et al.l[2000cj) . There is clearly a 
tight co lor-magnitude relation in t he cluster, fully ana- 
lyzed in Ivan Dokkum et alJ ljl998a|l . The entire catalog 
of spectroscopically confirmed cluster members within 
the HST mosaic is shown by the circles, and those se- 
lected for fundamental plane analysis are shown by the 
filled circles. As can be seen in the figure, the sample 
fairly represents the cluster population, minus the faint 
( an d very blue) objects . 

In lFisher et alJ ljl998() there are 108 members brighter 
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Fig. 2. — Overlayed on a model spectrum of a 6 G yr-old population JBruzual fc CharloHI20ri3) is the set of Lick/IDS bandpasses used in 
this paper IWorthev et al.lll994l:ITrager et aTHTOQ^ 



than R = 21 mag within the HST mosaic. Our high- 
resolution spectroscopic sample contains 52 of them. 
Three cluster galaxies fainter than the magnitude limit 
were added, bringing the tota l number of the samp le to 
55. The spectra were shown in lKelson et all ll2000bD . and 
the signal-to-noise ratios ranged from about 15 per A to 
90 per A. The 900/mm grating provided a typical reso- 
lution of ainst ~ 60 km/s. Details of the reductions to 
one-dim ensional spectra were also given in lKelson et al.l 

Below we describe the methods for measuring the ab- 
sorption line strengths and velocity dispersions. The al- 
gorithm for deriving the latter is necessary because we 
employ portions of it in determining the corrections for 
Doppler and instrumental broadening, and for estimat- 
ing the formal uncertainties in the line strengths. Fuller 
discussions of our measurements of the absorption line 
strengths, including the procedures to estimate formal 
errors, the procedures to correct the line strengths for the 
Doppler and instrumental broadenings, and the correc- 
tions for the di fferent metric aperture sizes are contained 
in .Kelson et al.i (,2006a.l . 

2.1. Absorption Line Strengths 

At this time the best understood and mod- 
eled diagnostics remain the spectra.1 indices of t he 
Lick/IDS svstem (iFaber et a l."1985: Burstcin ct al. 1984; 
Burstein. Faber. fc Gonzalea .1986; Gorgas et al. 1993 ; 
Worthev et al.lll994|) . defined as a set of indices mea- 
suring the strengths of Balmer, metal, and molecular fea- 
tures. We adopted the Tragcr ct al. ( 1998) definitions for 
the Lick/IDS indices, with the additional definitions of 
IWorthev fc Ottavianil l)1997t) for the indices of the higher 
order Balmer lines US and II7. Each index is comprised 
of a blue continuum bandpass, a central index bandpass, 
and a red continuum bandpass. 

The bandpasses of the blue indices are shown graphi- 
cally in Figure 121 Note that the stellar absorption lines 
covered by each index are heavily blended in galaxy spec- 



tra as a result of Doppler broadening and smoothed fur- 
ther by the line-spread function of the instrument. In 
order to compare line indices from one galaxy to an- 
other, or from data taken with one instrumental set-up 
to data taken with another, the indices must be corrected 
to identical levels of intrinsic Doppler and instrumental 
broadenings, and these corrections are discussed below. 

The absorption line strengths were measured using the 
following prescription: 

The straight line between the mean fluxes in the two 
continuum bandpasses is defined as the local continuum, 
C\, within the index bandpass. The mean flux within a 
bandpass is written as: 

Each index is then defined by the integral of the resid- 
ual flux above and/ or below the lo cal continuum within 
the index bandpass l|Worthevlll994D . For convenience we 
define 

The metal and Balmer line strengths have units of A and 
are commonly referred to as pseudo-equivalent widths, 
while molecular features are expressed as magnitudes: 

^ew=(A2-Ai)(l-/p) (3) 
^mag = -2.5l0g/p (4) 

The integrations are straightforward, with the contribu- 
tions from partial pixels computed using linear interpo- 
lation. 

2.2. Velocity Dispersions 

Many scaling relations of early-type galaxies use 
galaxy velo city dispersion, a, a, s the dependent vari- 
able (e.g. iFaber fc Jacksonlll976t iTerlevich et"al] I198H 
iDiorgovski fc Davislll987HDressler et alJll987j) . Rehable 
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measurement of velocity dispersion is also key to deriv- 
ing reliable line strengths as the velocity dispersions are 
required to compute the Doppler corrections to the line 
strengths. The procedures used for deriving a are also 
used in our method to estimate the formal errors in the 
line strengths. 

Our preferred method for measuring velocity dis- 
persions is t he dir ect-fitting method ( described by 
Kelson et all I2000bl though also see l|Rix fc Whitd 
19921) ). [t has several advantage s over te chnique s that 
operate in Fourier space (e.g . iTonrv fc Davia Il979t 
iFranx. Illingworth. fc Heckmad 119891) . Most impor- 
tantly, when measuring velocity dispersions of high- 
redshift galaxies, non-uniform sources of noise become 
important, and these affected pixels should not be ac- 
corded uniform weight. By fitting spectra directly, one 
measures the velocity dispersion by shifting and broad- 
ening a template sp ectrum until the result match e s each 
galaxy spectrum ijBurbidge. Burbidge. fc FishI 119611: 
IRix fc Whitd [199^ . Mathematically, one has a galaxy 
spectrum, G, and a template spectrum, T, and one 
searches for the line-of-sight velocity distribution i?, 
that is used to cony olve T such that is minimized: 
lIKelson et al.ll20n0bD : 

K 

= |{G- [PM(i?('T,t;)oT) + ^a,i/,]} xl^l'. (5) 

Following iKelson et al.l l)2000bf) we parameterize the 
linc-of-sight velocity distribution as a Gaussian, -B(cr, u), 
in which a is the velocity dispersion, or second moment 
of the velocity distribution and, v is the mean radial ve- 
locity, or first moment of the velocity distribution. The 
fit also includes two additional components: Pm-, a low- 
order multiplicative polynomial that incorporates the dif- 
ference between the instrumental response function in 
G and T, and ^ '^jUj is an additive continuum func- 
tion comprised of sin es and cosines up to wavenumber K 
ijKelson et al.ll2000b() . These sines and cosine continuum 
functions are equivalent to filtering low wavenumbcrs in 
the Fourier domain (see, e.g., Tonrv fc Davis 1979). 

The unique value of the direct-fitting method is con- 
tained in W , the w eighting spectrum. We opt, as in 
IKelson et all lj2000bD . to weight most pixels by the in- 
verse of the noise, estimated from photon statistics and 
amplifier read noise. Just as in Kelson et al. (2000b), 
there are exceptions: (a) regions contaminated by galac- 
tic emission lines are given zero weight; (b) strong Balmer 
absorption lines are also given zero weight; and (c) pixels 
contaminated by poorly subtracted night-sky emission 
lines are also given zero weight. Nearly all of the early- 
type galaxies discussed in this paper have no (or very 
weak) emission. The BCG, however, does have emission 
lines and therefore it is the only galaxy in this paper 
for which this "masking" is important. The fitting was 
performed between ^ 4250A and ^ 5050A, in the rest- 
frame, employed 5th-order for the multiplicative poly- 
nomial. Pa/, and si nes and cosines up to K = The 
results, as shown bv IKelson et alJ l)2000bl) . are sensitive 
to these choices only at the level of a few percent. 

The template spectra spanned a range of late-type 
stars, a s was also discussed previously in ' Kelson et al.l 
Il2000bl) . Ideally one might wish to employ spectra 
which represent more accurate models of the galaxy spec- 



tra. However, t he current state-of-the-art high-resolution 
model spectra l|Bruzual fc Gharloi]l2003l) have lower res- 
olution (^ SA) than our instrumental setup (~ 1.7A). 
Therefore these model spectra are not suitable for pro- 
viding accurate measurements of velocity dispersion in 
our lower-mass ga laxies. If one derives velocitv disper- 
sions using the Br uzual fc Chariot l)2003() SEDs, the one 
must add, in quadrature, the difference in resolution to 
the result. Doing so, we find < 1% systematic differ- 
ence bet ween those signias aii d the ones we obtained 
using the IKelson et alJ l)2000b|) stellar templates, with 
a standard deviation of 2%. Furthermore, there is no 
statistically significant trend with velocity dispersion for 
galaxies with tr > 100 km/s. 

2.3. Absorption Line Strength Errors 

Accurate estimates of the uncertainties in absorp- 
tion line strengths require one to determine the vari- 
ances with in each index bandpass (e.g., iGonzalez 199^ 
iGardiel et al. 1998). One can adopt the expected noise 
due to photon statistics and amplifier read noise, but 
sources of noise become increasingly localized for high 
redshift galaxy spectra, with the bright OH lines and, 
potentially, fringing dominating the errors in flux. Fur- 
thermore, the noise can become correlated in adjacent 
pixels after the data are rebinn ed onto a simple wave- 
length scale l|Gardiel et al.l ll998l 

We developed an approach that provides for a more 
robust determination of the line strengths and their er- 
rors, with less susceptibility to systematic errors. This 
approach requires neither the additional overhead of vari- 
ance images nor any assessment of correlated noise in ad- 
jacent pixels. Furthermore, our method does not rely on 
the assumption that the "true" noise is well-represented 
by the theoretical expectation from photon statistics and 
the electronics noise. Tests with the data and simula- 
tions indicate that our methodology produces accurate 
estimates of line strength errors. We discuss the method 
here. 

Deriving the random errors in absorption line strengths 
requires knowledge of the total variances in the contin- 
uum and index bandpasses. Given an observed galaxy 
spectrum G, the noise in a given pixel is simply TV = 
G — R, where R represents the underlying, noiseless spec- 
trum of the galaxy. While there is no a priori knowledge 
of i?. Equation fallows us to write 

K 

R^R = PM{BoT)+^ajHj (6) 
i=o 

Using the mod e l ste llar population spectra of 
IBnizua,] fc Gha.rlotJ we generated the best-fit 

model SED for each galaxy. The model spectra were 
used, instead of our template stars, because they bet- 
ter reproduce the Balmer, metal, and molecular line 
strengths in G. One should note, however, that 
in searching for the best-fit SED, no unique best-fit 
age/metallicity pair could be found for any given galaxy: 
young, metal-rich model SEDs fit as well as old, metal- 
poor model SEDs in all cases, with complete degeneracy. 
Because individual line strengths are subject to a similar 
degeneracy between age and metallicity (e.g., iWorthevI 
Il994|) . attempts to match the detailed spectra are sub- 
ject to the mean correlation between age and metallicity 
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because the computation of utilizes all of the spectral 
features simultaneously. 

The approximation to R provides an approximation to 
N KiN. 



K 



N = G~ [Pm{B oT) + J2 ajHj] (7) 

3=0 



With N the variance in the mean flux, Fp, is 
1 " ~ 



(8) 



where the summation is over the pixels within the given 
bandpass and n is the number of pixels within the band- 
pass. For pixels of constant AA, this reduces to the fa- 
miliar 

1 " ~ 
n '-^ 

i 

where the error in the mean flux is then 



\ 



1 " ~ 



(10) 



Note that partial pixels are included at the ends of the 
bandpass by linearly interpolating the fractional contri- 
bution of that pixel to the bandpass and adding those to 
the summation. For clarity we did not include them in 
the above equations. 

With accurate estimates of the errors in the mean 
flux of each bandpasses the total uncertainties in the 
line strengths can be trivially computed by propagat- 
in g the erro r s in a manner mathematically equivalent 
to iGonzale^ (|1993'). The key to our method is that 
the determination of the variances in each bandpass 
do not rely on simple photon statistics. In detail, the 
models may not provide perfect matches within the 
line strength bandpasses but these detailed departures 
from the models only become important when template- 
mismatch dominates the errors in the velocity disper- 
sion fitting. In the original reference for these spectra 
ijKelson et al.ll2000bj) . template-mismatch appears to be- 
come significant beyond S/N > 60 per A, but does not 
dominate the residuals between a galaxy spectrum and 
the model fit to it until one has significantly higher S/N 
ratios. Furthermore, by including functions in the fit 
for matching the continuum, broad mismatches are im- 
proved, such as for CN or C4668. 

The accuracy of our error estimates was confirmed in 
two ways. First, we performed Monte Carlo simulations, 
in which we degraded model spectra to a range of S/N 
ratios and remeasured the absorption line strengths. We 
found excellent agreement between our error estimates 
and the scatter in the measured line strengths about the 
input measurements. 

The second approach we took to verify the validity of 
our error estimates involved exploiting the fact that each 
spectrum had originally been observed with three expo- 
sures. For each exposure we measured the indices and 
their formal errors and created histograms of the distri- 
butions of the indices about their weighted means. These 



are shown in Figure|21 In the figure, Gaussians with stan- 
dard deviations of unity (the dashed lines) show the dis- 
tributions one expects if our formal error estimates are 
consistent with the differences between the three sepa- 
rate observations per galaxy. The solid lines show Gaus- 
sians whose standard deviations are given in each panel 
and were determined empirically from the observed dis- 
tribution of line strengths from the individual exposures. 
Excluding Mg&, the mean standard deviation of the his- 
tograms is 0.95, indicating that our errors, on average, 
are over-estimates by approximately 5%. Unfortunately, 
Mg6 is contaminated by telluric absorption at this red- 
shift and our corrections to it are uncertain. Thus for 
that index the variance between the measurements ap- 
pears to be systematically reduced, and, we infer, leads to 
an artificial r eduction in {x^) for its histogram (also see 
iKelson et al.ll200 6a. for examples of the spectra). Taken 
together, these tests indicate that our method for esti- 
mating the random errors in the absorption line strengths 
is accurate for spectra spanning the range of signal-to- 
noise ratios of these data. 

2.4. Correction for Instrumental and Doppler 
Broadening 

In spectra of galaxies, the blending of the absorp- 
tion lines significantly alters the continuum, and spectra 
with different intrinsic Doppler or instrumental broaden- 
ings will yield different pseudo-equivalent widths. Con- 
sequently, measurements of absorption line strengths 
must be corrected for both the internal motions of 
the stars and for the finite resolution of the spectro- 
graph. While the system of line strengths could have 
been defined by spectra with arbitrary levels of intrinsic 
Doppler or instrumental broadening, the Lick/IDS sys- 
tem of indices is referenced to stars (i.e., spectra with 
zero Doppler broadening) observed through the IDS, for 
which the instrument a l reso lution was characterized in 
iWorthev fc Ottavianil l|1997f l. It is important to note 
that this characterization of the IDS's instrumental res- 
olution has a very strong dependence on wavelength. 

The broadening corrections were computed by measur- 
ing an index, X, three times for each galaxy: {1) X\g, the 
index measured directly from the galaxy spectrum; (2) 
A^|s((t)oTj the index measured from the template spec- 
trum (or model SED; see previous section) broadened to 
the velocity dispersion of the galaxy and at the resolu- 
tion of the galaxy's spectrum; and (3) Xj^^^g, the index 
measured from the template spectrum (or model SED) 
with a velocity dispersion of zero and convolved to the 
resolu tion of the Lick/IDS system ijWorthev fc Ottavianil 
Il997|) . Thus, each galaxy's corrected line strength Xcorla 
can be recovered by 



cor I G 



X\tios ^ -'^|s(ct)oT 



(11) 



There are two advantages to defining the broadening cor- 
rections in this way. First, smoothing the galaxy spectra 
to low spectral resolution dramatically, and non-trivially 
alters the statistics because of non-uniform sources of 
noise. Second, our definition of the broadening cor- 
rections allows for sign reversals that can occur as the 
pseudo-continua change with large differences in Doppler 
and instrumental broadening. In these cases, such as for 
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Fig. 3. — Histograms showing the variance in the measurements of the indices from the individual exposures. Each galaxy's spectrum 
was observed in three exposures and a comparison between the line strengths in the individual exposures provides an independent test 
of our estimates for the formal errors. The dashed curves show Gaussians with standard deviations of unity to illustrate the expected 
distribution of the data. The solid curves show Gaussians with standard deviations determined empirically and listed within each panel. 
Excluding Mg6, the mean standard deviation is 0.95, indicating that the formal errors in our line strengths have been properly estimated. 
The Mgb bandpass suffers from telluric absorption and our correction to it is uncertain. As a result the residual telluric absorption reduces 
the variance between the individual observations. 



the h igh-order Balm e r hnes, mult i plicat ive corrections 
(e.g., IGonzalezllT993t iTrager et all [1991 can be prob- 
lematic. 

Because broadening corrections can be sensitive to 
the adopted temp late SED, we adopt the best-fit 
iBruzual fc Charloti pOO^ SED for T. By doing so the 
systematic error in the broadening corrections is mini- 
mi zed. However, the ages an d metallicities of the best- 
fit IBruzual fc Charloti i)2003{) SEDs are sensitive to the 
wavelength range of the fit. Therefore we derive broad- 
ening corrections using parameters derived from fitting 
over 4100A < A < 5000A, 4200A < A < 4800A, 
4070A < A < 4800A, and 4070A < A < 4500A, in 
the restframe. We adopt the mean broadening correc- 
tion from these four fits, and we use the scatter in their 
broadening corrections as estimates of the random errors 
in the corrections. 



Fortunately, the large uncertainties in the SED ages 
(~ 0.15 dex) and metallicities (~ 0.3 dex) are strongly 
(anti)correlated such that the uncertainties in the broad- 
ening correction are small. For all of the indices, the 
errors in the broadening corrections are estimated to be 
at a typical level of ±10% of the correction itself; i.e., 
quite small. More importantly, the uncertainties in the 
corrections are typically < 3% of the dynamic ranges of 
the line strengths, as estimated using the 68% widths of 
the distributions of the indices. 

Our broadening corrections are plotted against galaxy 
velocity dispersion in Figure 01 At velocity dispersions 
(7 < 150 km/s, there is an increase in the scatter in the 
broadening corrections for several indices. Even though 
these increases in the scatter arise because of an increase 
in t he scatter of the best-fit a ges and/or metallicities of 
the IBruzual fc Chariot 1)20031) SEDs, the increased scat- 
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Fig. 4. — Broadening corrections to the blue Lick indices for the E/SOs in CL1358+62 as a function of velocity dispersion. Symbols are 
as i n Fipj. pi The increased sc atter at low velocity dispersions arises because of an increase in the scatter of the ages and/or metallicities of 
the'Bruzual fc C harloU 1200? ) SEDs used to derive the corrections for each of the galaxies. Much of this increased scatter is not statistically 
significant, though the variations in the intrinsic strengths of the features, i.e. variations in the underlying SEDs, are important for indices 
such as H7^. 



ter in the corrections is typically not statistically signif- 
icant, given the uncertainties in the corrections them- 
selves. A small number of galaxies fall off the primary 
loci of points for the narrow Balmer indices HJi? and 
H7i?. These galaxies, below the magnitude limit of the 
sample, do have larger values of these indices, though the 
wider H(5^ and indices, which will be used in our 
analysis below, do not show such large departures. 

2.5. Correction for Aperture Size 

In order to accurately compare samples of galaxies at 
different redshifts, one must be correct for the effects of 
observing the galaxies with different aperture sizes. Be- 
cause this sample will be the reference datase t in our sur- 
vey of distant chisters [Kelson et alJ l)2006bD . we correct 
the observed line strengths to an effective aperture con- 
sistent with the aperture used in the most d istant cluster 
in our sample, MS1054-03, at z=0.83. In Kelso n et al.l 
we detail our formalism for estimating the re- 
lationship between the internal line strength gradients 



in distant galaxies and the dependence of line strength 
measurements on aperture size. 

Using estimates of the line strength gradients for the 
CL1358-I-62 sample galaxies, and the relationship be- 
tween gr adient, 0, and apert ure correction, AX, de- 
rived in 'Kels on et alJ l)2006afl . we have corrected the 
line strengths of the CL1358-I-62 galaxies to an aperture 
equivalent to a diameter Z?iof — 10.1 kpc at the dis- 
tance of MS1054-03 {Dap = 1''35, using, again, HQ = 72 
km/s/Mpc, = 0.27, and 17a = 0.73). The correc- 
tions are defined to be AX = (3 x fiDgn/Drei ) , wher e 
we derived f{Dap/D,ei) = 0.075 in lKelson et all l|2006al) . 
The corrections for aperture size are small because the 
extraction aperture for the CL1358-I-62 galaxies was cho- 
sen to minimize the corrections. As a result, the errors 
in correcting the aperture to an equivalent aperture at 
the distance of MS1054-03 are smaU (< 10% of the cor- 
rection). 

The uncorrected line strengths for the full sample of 
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galaxies in CL1358+ 62, are given in several tables of 
iKelson et al.l l)2006aD . The line strengths of the ellipti- 
cal and lenticular galaxies, corrected for both broadening 
and aperture size, are reproduced here in Tables and 
IA2I Though these are the data to be used in the discus- 
sions of the stellar populations below, these corrections 
do not affect the conclusions drawn in this paper. 

2.6. The Line Strength-Line Width Relations 

The correlations between absorption line strengths and 
velocity dispersion has become a diagnostic for d i rectly 
probing galaxy evolution (e.g. I Zicglcr & Be nded 119971 
IKelson et a'l]l200lHJorge nsen et al.i2005^1 . The strengths 
of the metal, molecular, and Balmer absorption features 
all correlate with velocity dispersion, and the full multi- 
dimensional locus of early-type galaxies in the space 
of velocity dispersion and line strengths has implica- 
tions for scenarios of their formation, e.g., such as the 
timescale for the build-up of stellar mass. For example, 
the (H74-|-H < 5/i)-g relation was studied to 2; = 0.83 by 
IKelson et al.l l)2001|) in an effort to directly trace the evo- 
lution of the high-order Balmer lines, at fixed velocity 
dispersion, with redshift. 

Before proceeding to detailed modeling of the indices, 
we show the observed line strength-line width relations 
for the early-type galaxies in CL1358-)-62 in Figure 
Note that the E and SO galaxies follow statistically iden- 
tical line strength-line width relations. This fact has 
consequences for the uniformity of their formation and 
enrichment histories, and we discuss this further below 
using simple models. 

3. FITTING FOR THE RELATIVE STELLAR POPULATION 
PARAMETERS 

The model s pub lished in 

Thomas. Maraston . & B endeij l)2003aH and 

Thomas. Maraston. & Korn (2004) contain predic- 
tions for the Lick/IDS indices at locations in sparsely 
sampled "grid" of the six parameters: \ogt, [Z/H], [a/Fe] 
(mean a enhancement), [a/N] (nitrogen enhancement 
or depletion on top of any enhancement of N as an a 
element), [a/C] (same as [a/N] but for carbon), and 
[a/Ca] (for calcium). For clarity, when nitrogen is 
enhanced with respect to the ensemble of a elements, 
[a/N] decreases. As a reminder, most published results 
on a-enhancement stem from the measurement of indices 
sensitive to Mg. However, oxygen is the most abundant 
a element and will likely dominate the a sensitivity 
in the blue owing to the equilibrium between carbon, 
nitrogen, and oxygen. 

The dependence of these spectral indices on the con- 
trolling parameters is non-linear and the inversion of ob- 
servables to find the luminosity-weighted mean proper- 
ties of stellar populations is non-trivial. In this section 
we detail the steps by which we determine the stellar 
population model parameters that best fit the observed 
line strengths. 

Each of our galaxies in CL1358-I-62 has between 
eight and ten indices measured from its spectrum. For 
every galaxy, we wish to use all of the available data to 
constrain its star-formation and nuclcosynthctic history, 
as described by the "simple" models. In principle 
our fit of a simple stellar population (SSP) to a given 
galaxy has up to six unknowns, and of order ten 



observables, suggesting that the process of inverting the 
observables should be straightforward. However, each 
Lick index measures the combined strength of many 
blended features in the spectrum of the underlying, 
unbroadened stellar population. These features are all 
sensitive to the ages, mean meta llicities, and d e tailed 
abundance ratios of the st ars iBurstein et alJ 119841: 
Worthev et al.llT99l lWorthevlll994l: iTraeer et al."2000bl: 
Cardiel et alJ 120031 iThomas. Maraston~fc BendeJ 
'2003at iThomas. Maraston. fc KornI |2004|) . To help 
illustrate these degeneracies, Table IA3I lists the first 
partial derivatives of the model predictions (at a specific 
reference point, discussed below). These derivatives 
indicate that the process of fitting models to observed 
indices will lead to significant covariances between the 
fitted parameters. The derivatives given in the table 
are only valid at a single location in the six dimensional 
space of SSP parameters and the full dependencies of the 
model predictions on the SSP parameters is nonlinear. 
Below we describe how the nonlinearity of the models 
has been treated. 

At this time we note that some derivatives are not 
listed. These are explicitly assumed to be zero, and 
we elaborate on the implications of this below, in the 
context of our method for fitting differential stellar 
population parameters. The unknown derivatives have 
occurre d because either IThomas. Maraston. fc Bended 
()2003a|) specifically indicated that that particular sen- 
sitivity was negligible, or the models do not yet include 
the sensitivities of these indices to those abundance ra- 
tios (foipicco^_&_^elIJ_995^^ & Koi3 
'2004 0Korn. Maraston. fc Thomasli2605l(l . The next sec- 
tion discusses our treatment of these shortcomings. 

In principle the wavelength coverage in our spectra al- 
lows us to utilize all of the blue Lick indices from USa 
through Fe5015 in the fitting of stellar population pa- 
rameters for most of the galaxies. Mg5 was observed for 
many galaxies but we do not use it because of the un- 
certainties introduced by the strong telluric absorption. 
The spectra of some galaxies also missed portions of the 
US and CNi continua or index passbands. In theory one 
could derive ages and metallicities for all of the galaxies, 
despite the absence of a few random indices from some 
galaxies. However, the dependencies of the parameters 
on the observables would not be uniform across the sam- 
ple and the resulting ages and abundances would not be 
homogeneous. 

In an effort to construct an homogeneous set of stellar 
population parameters, we opt to fit only those galax- 
ies for which there exist measurements of ^ISa, ^IAi 
CN2, Ca4227, G4300, Fe4383, Fe4531, and C4668 (H^ 
has been excluded from the fitting because it can be af- 
fected by emission). Out of the original sample of 55 
ga laxies, 33 of which were classified as E, E/SO, or SO 
bv'Fabricant^ Franx^ fc v an DokkumI l|2000(l in the orig- 
inal Kelso n et alJ l)2000b|) sample, we are left with an 
homogeneous sample of 19 early-type galaxies for the de- 
tailed study of stellar populations. Our primary analysis 
will rely on this homogeneous sample. To help illustrate 
whether the fuller sample has comparable uniformity to 
the smaller, homogeneous sample, the remaining galaxies 
are included in some of the figures. As can be seen in Fig- 
ure 13 these additional objects obey the same relations 
as the smaller subset, with no significant outliers. 
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Fig. 5.— The line strength-line width relations for E/SOs in the CL1358+62 sample of Kelson et al.' f 2000c'). Ellipticals are shown as 
filled red circles and E/SO and SO galaxies as open blue circles. Note that the measurements of Mgfe, not used in the fitting for the stellar 
population parameters, suffer from additional uncertainties arising from poor correction to telluric absorption. The small symbols represent 
those galaxies that are missing any of the eight required indices for the fitting, and as such belong to a larger inhomogeneous sample. The 
larger symbols show those galaxies which have all eight indices to be used in the stellar population fitting. 



3.1. Eliminating the Zero-point 

Before fitting the indices to find the relative ages and 
abundances, one more issue requires resolution: system- 
atic offsets between the system of measured indices and 
the mo del system, defined to be equivalent to the Lick 
system ffiurstein et"allll984t iTrager et al.lll998j) . There 
are two distinct issues here: (1) actual differences be- 
tween the observations that form the basis of the sys- 
tem of model indices and the indices measured from 
our spectra; and (2) unknown dependencies of indices 
on any of the abundance ratios. We have attempted to 
correct for differences between the model and observed 
systems of indices using the procedures described ear- 
lier, but significant systematic offsets between the model 
predictions and the data may remain. Such systematic 
differences between the systems of the models and our 
data severely hamper our ability to infer the absolute 
ages and chemical make-up of the stellar populations in 
our sample. Separate from differences in the systems 
of the modeled and observed indices, intrinsic calibra- 
tion uncertainties exist for the model indices themselves. 
For example, the higher ord er Balmer lines defined by 
iWorthev fc Ottavianil 1)19971) . may have calibration, or 
zero-point uncertainties of up to ~lA llTragei<l2002t) . 

The other systematic effect, that of ignored, or even 
unknown dependencies of the line strength predictions on 



any of the various abundance ratios, is also very impor- 
tant. Recently Thomas. Maraston. & Korn ( 2004) have 
determined the sensitivities of the high-order Balmer 
lines to [a/Fe], and those authors argued that past 
discrepancies between ages derived from the different 
Balmer lines can be attributed to the previously uncal- 
ibrated dependencies of those features on [a/Fe]. How- 
ever, other elemental dependencies, such as a dependence 
of HjA on nitrogen or carbon , are also likely to be im- 
portant ijSchiavon et al.ll2'002() . 

We eliminate these two sources of systematic error by 
zero-pointing the models to our data. Doing so ensures 
that every index yields the same age, metallicity, etc. 
This serves to eliminate any zero-point uncertainties in 
the indices, arising from calibration errors or from un- 
known zeroth-order dependencies of the line strengths 
on all of the abundance parameters. The chief penalty 
incurred by recalibrating the models is the loss of the 
zero-points for logi, [Z/H], [a/Fe], etc.. More specif- 
ically, the model's stellar population parameters are no 
longer absolute but are relative. Fortunately we are most 
interested in deriving accurate relative measures of age 
and abundance, and the absolute values are simply not 
important for this paper or for our survey of the proper- 
ties of cluster galaxies as a function of look-back time. 

The recalibration is performed by adopting a set of 
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stellar population parameters for a reference point in 
the space of observables. This reference point is defined 
by the mean massive early-type galaxy in CL1358-I-62, 
computed, for each index, using the means of the 
line strengths for the most massive E/SOs in the clus- 
ter. We adopt a mean age of 7 Gyr (i.e., Zf — 
2.4: iKelson et alJ 120011 for the reference point, based 
largely in the results from previously published stud- 
ies of the evolution of E ^SO ga laxies. (Kclson ^et al .l 
1997; 'van Dokkum et al.' 1998b; Kelson et al. 2000c| 
" n Dokkum fc Stanford 2003; Holdcn et al. 2004, 200^ 
and others). We adopt abimdanccs of [Z/H] = 0.3; 
[a/Fe] = 0.2; and [a/N] = [a/C] = [a/Ca] = because 
work on nearby ellipticals has shown that the most mas- 
sive of such galaxies typically have super-solar metallicity 
and g-enhancement s between 0.2 < [a/Fe] < 0.3 (e.g., 
iTrager et al]l2000b|) . Within the range of values defined 
by previous work on nearby galaxies, our key results are 
not sensitive to changes in the adoption of these specific 
parameters. In Appendix we show that our key re- 
sults are insensitive to uncertainties in the adopted SSP 
parameters of the recalibration. 

The recalibration of the models is performed in the fol- 
lowing fashion: (1) model line strengths are generated for 
the reference SSP given above; (2) mean line strengths 
are computed using the 9 E/SO galaxies with a > 200 
km/s; and (3) arithmetic differences between these mean 
line strengths and the reference model line strengths are 
then used as zero-point corrections to the models when 
the fitting is performed. By doing so, we ensure that the 
data will always fall within the bounds of the model grid, 
and the location of the minimum for a given galaxy 
will be easily accessible in a downhill search beginning at 
the reference point. As has been pointed out by other au- 
thors, the largest systematic error in the location of the 
minimum arises from uncertainties in the zero-points 
of the models (vs one's observations). By eliminating 
this major, non-Gaussian source of uncertainty, the fit 
for the relative population parameters can now be per- 
formed by minimizing x^, with the expectation that now 
the covariance matrices can be used to estimate formal 
uncertainties. 

This last point is important: even if our assumed zero- 
point stellar population parameters are wrong (and they 
probably are) , the derivatives in the 6-dimensional model 
grid will be correct to first order, and thus the relative 
offsets Alogi, A [Z/H], will be correct to second or- 
der. If we had assumed that the models and our data 
are on the same system, then large systematic errors in 
the calibration would have translated into large uncer- 
tainties in the absolute values of the stellar population 
parameters. Furthermore, because different indices have 
different (or unknown) dependencies on several of the 
abundance ratios, these systematic uncertainties would 
have led to large errors in the ages and abundances. 

Essentially such systematic uncertainties in the zero- 
points of our measurements result in different absolute 
ages or abundances for each index. By adjusting the 
zero-points of the six-dimensional grid for a given index, 
one eliminates the zeroth order term in the dependence 
of that index on the abundance ratios. For indices whose 
dependencies on abundance ratios have been modeled, we 
know that the derivatives vary slowly, and that the mixed 
second derivatives are small. Taken together, the errors 



in the relative stellar population parameters essentially 
have only second-order dependencies on the abundance 
ratios. 

If the total systematic error between the calibration of 
the observations and models had been ignored, we would 
not be able to correctly estimate formal uncertainties, 
even for relative measures of age and abundance. Bo- 
gus correlations between the stellar population param- 
eters would have appeared, and underlying correlations 
between parameters might have been masked. These is- 
sues are made more explicit in Appendix^ All of these 
issues were also explicitly verified when we incorporated 
the upgraded models, in which the [a/Fe]-dependence of 
the H7 and US indices were included, and our results on 
the differential properties of the stellar populations did 
not change in any statistically significant way. 

The model line strengths at the reference point are 
given as Xmodci in Table IX4l The mean (adopted) line 
strengths for the reference point are given as ^(CLisss) • 
The difference, or zero-point corrections, are given by 
Ax- Again, the Appendix 1X1 discusses the impact of 
ignoring this zero-point recalibration. 

3.2. Minimization 

The derivation of the stellar population parameters for 
a given galaxy is performed by searching through the 
six dimensional space of parameters to find the set that 
minimizes = I] [(^modci + Ax - X)laxf (also see 
IProctor et a,ni200l . where X G { H5^, CN2, Ca4227A, 
G4300A, H7A, Fe4383A, Fe453lA, and, C4668A}, ax 
is the formal error in AT, Xmodei is the model predic- 
tion for index A", and Ax is the zero-point correction 
for the index, as described above. In order to linearize 
the minimization of x^, smooth approximations to the 
model grids are required, and to this end the discrete 
model grids were replaced with multivariate B-splines. 
The model indices are fit by bicubic B-splines in both 
\ogt and [Z/H], where the knot coefhcients are them- 
selves represented as low-order polynomial functions of 
[a/Fe], [a/N], [a/C], and [a/Ca]. Constructed in this 
way, the original multidimensional grids are reproduced 
exactly, and the interpolating function has well-behaved 
derivatives everywhere the bivariate B-spline is defined: 
< logt < 1.3 (G yr) and -2 < [Z/H] < 1. The 
iThomas. Maraston. fc Bender (2003a) models are de- 
fined between —2.25 < [Z/H] < 0.67. Fortunately none 
of the sample discussed here have fitted [Z/H] values out- 
side of this range. However, any galaxies that would have 
been fit using very large metallicities, e.g., [Z/H]> 0.67, 
would have properties extrapolated from the grid. 

The x^ minimization is performed using 

a modified Levenberg-Mar quardt algorithm 

ijGarbow. Hillstrom. fc Mord il980() . For single stellar 
populations, the mapping of parameters to observables 
is unique. Therefore, we begin each search at the 
reference point, at old ages, solar metallicity, and solar 
abundance ratios. From that portion of the grid, the 
best-fit stellar population parameters are accessible by 
any efficient gradient search algorithm. 

The result of a fit is a set of stellar population parame- 
ters defined relative to the mean of massive E/SO galax- 
ies in the cluster. These differential SSP parameters are 
referred to by Alogt, A [Z/H], A [a/Fe], .... Note that 
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we are fitting all of the available data, while others (e.g. 
iCardiel et al.ll2003^ have suggested that one can reduce 
the number of observables to include only those that are 
most sensitive to the desired stellar population parame- 
ters. Our methodology does not require such an analysis. 
Because we have eliminated the systematic uncertainty 
in the zero-points, and have reliable estimates of the for- 
mal errors in the line strengths, the Jacobian correctly 
accounts for the sensitivity of the parameters to all of 
the measurements (normalized by their uncertainties). 

4. THE RELATIVE AGES AND CHEMICAL ABUNDANCES 

The properties of the stellar populations of the 19 
early-type galaxies are found by fitting to only eight of 
the available observables. And while the models do allow 
one to fit for six parameters, it would be foolhardy to do 
so, given the uncertainties in our measurements and the 
degeneracies in the models. Therefore we employ a se- 
ries of constrained models. The results of the fitting will 
be discussed within the context of finding the minimum 
number of parameters required to define the family of 
early- type galaxies. 

4.1. Models with Variable Ages 

Here we explore several sets of models in which the 
relative ages are fit, and in which the parameters con- 
trolling the relative abundance ratios are optionally kept 
frozen. These models will be used to explore to what 
extent the SSP ages vary among the early-type galax- 
ies in this sample. Figure El^a) shows results when 
one fits only for Alogt and A[Z/H], and keeps the 
other four parameters fixed. The filled red circles de- 
note the ellipticals; the open blue circles the E/SO and 
SO galaxies. The morphologies, a gain, are taken from 
iFabricant. Franx. fc van DokkunJ l)2000|) . Because the 
models have been recalibrated to the data, each index 
should yield the same ages and abundances. However, 
the different indices have different sensitivities so data 
missing for particularly indices can bias the results for 
those galaxies. Therefore, only those galaxies in the ho- 
mogeneous sample, with all eight observables, are shown 
and all discussions of correlations and statistical varia- 
tions will refer only to these galaxies. 

In each of the figures for this section, we plot the best- 
fit relative stellar population parameters against galaxy 
velocity dispersion, and also include a panel showing the 
reduced of the fit for each galaxy. In the plots of rel- 
ative age, or metallicity, vs. velocity dispersion, we show 
the approximate location of magnitude selection cut us- 
ing three gray tones. The darkest gray region shows the 
region of velocity dispersions for galaxies below the mean 
dispersion at i? = 21 mag, using the Faber- Jackson re- 
lation of the full sample. The -1-1- and +2-a scatter of 
the i?- velocity dispersion correlation is indicated by the 
lighter gray strips. In these regions the galaxy selec- 
tion is no longer random and the sample becomes biased. 
Therefore, we can not include those galaxies in any cal- 
culation of correlations of the parameters with redshift. 

In each figure we show the effects of the selection us- 
ing the thin and thick solid lines, derived by fitting cor- 
relations using those galaxies with a > 134, 158 km/s, 
respectively. These cuts in velocity dispersion represent 
the mean dispersion, and the mean dispersion plus one 
standard deviation, at i? = 21 mag. In principle, 50% 



and 32% of the underlying sample of galaxies at the mag- 
nitude limit of i? = 21 mag are lost. Those galaxies at 
the limit which are over-luminous for their velocity dis- 
persion, either because they are younger or metal poor, 
can make it into the sample and bias the observed pop- 
ulation trends. In fitting for correlations, we opt not to 
weight using the uncertainties because the measurement 
errors are significantly correlated with a. Uncertainties 
in the fitted slopes were computed using the "bootstrap" 
method (with 5000 random samples). We discuss any 
specific correlations in more detail below. 

For the models in which only Alogt and A[Z/H] vary, 
the mean per degree of freedom for the two-parameter 
models was (x^) = 1.86 ^ 1, indicating that such mod- 
els, in which only these two parameters vary, are not 
complete descriptors of the stellar populations. How- 
ever, despite the poor quality of the fitting, the [Z/H]-ct 
relation is apparent. Discussions of the significance of 
any correlations with a are deferred until below, when 
fits to the data with lower (x^) are derived. 

In fitting the line strengths we find that the two pa- 
rameters Alogi and A[Z/H] are insufficient for describ- 
ing the observations. Therefore, we proceed to add ad- 
ditional parameters, one-by-one. Because of the previ- 
ously published work incorporating a/Fe as a key pa- 
rameter in the stellar populations of earl y-type galaxies 
(e.g.,D 0rgcnsc l3ll999l:lTrager et a l. 2000b: iProctor et alJ 
l2004ariThomas et al.ll2005D . we now allow A [a/Fe] to 
vary in the fitting. In Fig. Efb), these results are shown, 
again plotted against velocity dispersion. For purposes 
of illustration, the gray regions of exclusion are plotted 
assuming that the M/L ratios have zero dependence on 
the abundance ratios. With (x^) = 1.92 it appears that 
a-enhancement, by itself, does not serve to improve the 
fits to the blue Lick indices. Furthermore, it suggests 
that a/Fe is not an important parameter for defining 
the sequence of E/SOs in CL1358-f62. The ages appear 
to be uniform and the A[Z/H]-log(T correlation is only 
slightly fiatter than in Fig. Efa). These models indicate 
that a/Fe is a constant in these galaxies, and not a strong 
function of velocity dispersion. 

These models shown in Fig. Elb) are directly com- 
parable to many previously p ublished studies o f early- 
type galaxies re.g.lKuntschneillT99 8: Tra ger et alJl2000bl: 
iThomas et al.ll2005HZiegler et al.ll2005l|) . A linear least- 
squares fit to the relative a-enhancement as a function 
of log cr gives a slope of —0.01 ± 0.18 dex/dex (the thick 
line in th e figure), com pared to the [a/Fe] oc 0.33 logo- 
found bv iTrager et alJ (2000b) or [a/Fe] oc 0.28 logo- 
found by ITho mas et all (,2005.1 . While we do not find 
a statistically significant correlation of A [a/Fe] in our 
3-parameter models, the difference between our result 
and theirs is discrepant at the 2-o- level (even without 
knowing what the uncertainties in their slope estimates 
are, and with our comparatively small sample). Further 
discussion comparing our constraints on [a/Fe] with pre- 
vious results is reserved for i )5.1l Nevertheless, adding 
A [a/Fe] to the fitting of the blue indices did not reduce 
(x^); surely other abundance ratios a r e more important. 

The IThomas. Maraston. fc Bended l)2003a(l models al- 
low three other abundance ratios to be varied. In pre- 
viously published analyses of line strengths, there has 
been evidence for variable nitrog en abundance ratios 
among old stellar populations, and l)Burstein et al.llT984l: 
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Fig. 6. — Stellar population parameters as a function of velocity dispersion for ellipticals (filled red circles) and E/SO and SO galaxies 
(open blue circles) for three different classes of models; (a) where only relative age and metallicity are fit (A log t-A[Z/H]); (b) where 
relative age, metallicity ,and o-enhancement are fit (A log t-A[Z/H]-A[o/Fe]); and (c) where relative age, metallicity, o-enhancement, and 
nitrogen-enhancement are fit (A log t-A[Z/H]-A[Q/Fe]-A[Q/N]). Only those galaxies with all eight observables are shown. The left-most 
panels show the reduced fo'' each galaxy, also plotted against galaxy velocity dispersion. The mean reduced values are also given 
in these panels. There are several key features to note: (1) There appears no statistically significant correlation of age with velocity 
dispersion; (2) there is a clear A[Z/H]-log(T correlation; and (3) our correlations of A[a/Fe] with logcr are marginally significant and they 
even appear to be anti-correlated when nitrogen is allowed to vary; and (4) The addition of A[a/N] as a fitted parameter significantly 
improves the fits to the observed line strengths. The darkest gray region shows the location of the R = 21 mag selection cut, using the 
i?-band Faber-Jackson relation in the cluster. The lighter shades of gray show the and +2-a scatter about the Faber-Jackson relation. 
In plotting the abundance ratios, we assume the M/L ratios have zero dependency on the specific abundance ratio in order to show the 
selection. Within these gray regions, the selection effects become important and the sample becomes non-random. As a result, we do 
not include those galaxies within these biased regions when fitting for population trends with tr. The thin and thick solid lines show the 
correlations obtained when fitting using galaxies with velocity dispersions greater than the mean dispersion at R = 21 mag (134 km/s), 
and the mean plus one standard deviation (158 km/s), respectively. The abundance trends are insensitive to the selection biases, while 
apparent correlations of age with velocity dispersion are biased strongly by those galaxies at the limit of the sample. As a result, there is 
no significant evidence for a correlation of stellar population age with a. 



IBrilev. Cohen, fc Stetsol]l200a l2004aD among ear ly-type 
galaxies specifically fe.g. iCarretero et al.l 1200^ . We 

therefore investigated the impact of the A[a/N] relative 
abundance ratios, and the results are shown in Figure 
EJc). First, note the significantly lower mean per de- 
gree of freedom, (x^) = 1.16, indicating that nitrogen 
is a key factor in fitting the blue Lick indices (in par- 
ticular the cyanogen index). Second, there is a strong 
correlation of the nitrogen abundance ratio with galaxy 
velocity dispersion. The origins of this correlation will 
be explored in a later section. 

In the models shown in Figure El the slope of the cor- 
relation of relative age with velocity dispersion is very 
sensitive to selection biases. Excluding galaxies below 
a < 158 km/s, one obtains Alogt cx (0.06 ± 0.19) loga. 
If one cuts the sample at the mean dispersion at i? = 21 
mag (ct < 134 km/s), one would expect that galaxies 
which are "young" for their mass/dispersion would make 
it into the sample, while older galaxies drop out, thus 
producing a bias of younger ages at lower velocity dis- 
persion. In fact, this bias appears to be present, as one 
finds a slope of Alogt cx (0.38 ± 0.18) logcr. This bias 
is not minimized by selecting only those galaxies with 
velocity dispersions greater than the mean dispersion at 
the magnitude limit of the sample, but by selecting only 



those galaxies with dispersions greater than the mean 
plus at least one standard deviation of the relation that 
transforms the magnitude limit to a cut along the ab- 
scissa. The larger the scatter in the transformation of 
the magnitude cut to a cut in the dependent variable, 
the greater the bias in any correlation is. 

Because of the above arguments, we conclude there is 
no significant evidence for a correlation of stellar pop- 
ulation age with a, down to the limit of our sample. 
Furthermore, we find no significant difference between 
the ages of the SO and elliptical galaxies down to our 
magnitude limit. Note that this result is consistent with 
what was found in ou r analysis of the fundamental plane 
(IKelson et al.ll2000c]) . 

The best fits to the data in which ages are allowed 
to vary are shown in Fig. Efc). These data indicate 
that the stellar populations of the early-type galaxies 
not only have remarkably uniform ages at fixed veloc- 
ity dispersion, but that they also follow a tight A[Z/H]- 
logcr correlation. The resulting distributions of Alogt, 
A[Z/H], A[a/Fe], and A[a/N] with velocity dispersion 
can be summarized by the following: 
• The scatter in relative ages, above the limit of the 
sample, is 0.06 dex. The scatter expected from the for- 
mal errors is 0.09 dex, suggesting that the formal errors 
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Fig. 7. — Same as in Figure 151 but where ages are assumed to be constant, (a) parameters from fitting only for A[Z/H], A[ci/Fe], and 
A[a/N]; (b) parameters from fitting only for A[Z/H], A[a/Fe], A[a/N], and A[a/C]; and (c) parameters from fitting only for A[Z/If], 
A[a/Fe], A[a/N], and A[a/Ca]. No improvements are found by allowing the Ca and C abundance ratios to vary. We find no significant 
correlation of these abundance ratios with veloc ity d ispersion. The ramifications of the correlations of the nitrogen and carbon abundance 
ratios with velocity dispersion arc discussed in il5.3l 



on the relative ages have been overestimated by a third. 
The observed scatter is consistent with th e low scatter 
in mass -to-light ratios, at fixed a, found in 'Kelso n et al.l 
()2000d) . There is no significant correlation of relative age 
with velocity dispersion, down to the limit of the sample, 
consistent with the lack of evolution of the fund amental 
plane slope to this redshift ijKelson et alJl2000^ . 

• No difference in the relative ages and abundance 
patterns of Es and SOs is seen in CL13 58+62, con- 
trary to what was seen in Fornax by iKuntschnerl 
1)19981) . This is presumably because the Es and SOs 
in our sample span simil ar ranges o f veloci ty disper- 
sions, while the SOs in IKuntschnerl l)1998|) all have 
(T < 100 km/s, and, as such, belong to a larger, gen- 
erally inhomogeneous set of early-ty pe galaxie s (e.g. 
Jorgensen. Franx. fc Kisergaardl Il995t ICaldwell fc Rose 
1997 1 iKelson et al.ll2000d iCaldwell. Rose, fc Concannon 

mm . 

• A steep, tight correlation of A[Z/H] with velocity dis- 
persion is found, with a slope of 0.86 ± 0.17 dex/dex. 
This is about twice as steep as that seen in other sam- 
ples of early-type galaxies (e.g. iTraeer et al.ll2006bj) and 
we investigate the implications for the color-magnitude 
relation below. The scatter in A[Z/H], at fixed a, is 0.06 
dex, where the scatter expected from the formal errors 
is 0.05 dex. 

• About 30% larger scatter is seen in A[a/Fe] than 
the 0.04 dex expected from the formal errors. There 
appears to be a mild anti-correlation of A[a/Fe] with 
velocity dispersion, with a slope of —0.44 ± 0.25 
dex/dex. Using those galaxies with a > 134 km/s, 
one obtains —0.28 ± 0.15 dex/dex, indicating poten- 
tial sensitivity to the selection. The inclusion of 
the nitrogen abundance ratio in the fitting does have 
an impact on the best- fit values for A[a/Fe], with 



many of t he b lue indices providing leverage on [g/Fe] 
(see Table IX3l and iThomas. Maraston. fc Benderll2003at 
iTh omas. Maraston. fc KornI \2004 ~ 

Without variable 

nitrogen abundance ratios, one obtains A[a/Fe] oc 
(—0.01 ±0.18) logCT. Such a-enhancement patterns differ 
from previously published correlations of A[a/Fe] with 
velocit y dispersion (see, e.g.lGonzalezlll993 : Traacr et alJ 
l2000bl:lMehiert et al.ll200a iThomas et al.ll2005^1 . though 
earlier work relied on redder Lick indices (e.g. Mg6), em- 
ployed models in which the nitrogen was not free to vary. 
We discuss the implications of these A[Q!/Fe] estimates 
on the Mg-fT relation below. 

• A steep anti-correlation of A[a/N] with velocity dis- 
persion is found, with a slope of — 1.36 ± 0.32 dex/dex. 
Using those galaxies with cr > 134 km/s, one obtains 
— 1.17 ± 0.22 dex/dex. Given the results for a/Fe, we 
can restate the correlation of nitrogen enhancement with 
velocity dispersion as [N/Fe] oc (0.94 ± 0.24) log cr, or 
[N/H] oc (1.78 ± 0.25) loga. The imphcations of the ap- 
parent nitrogen-(T correlation will be discussed further 
below, though at a basic level the nitrogen abundance 
ratio appears to play an important role in dictating the 
behavior of many blue Lick indices, especially given its 
impact on the fits for [a/Fe]. The scatter in A[a/N], at 
fixed (7, is 0.08 dex where the scatter expected from the 
errors is 0.06 dex. 

4.2. Models with Variable Carbon and Calcium 

The results of the fits in the previous section indicate 
that the ages of the early-types show no statistical vari- 
ation down to the limit of our selection. Therefore we 
now fit models in which we forced Alogi = 0. Fixing 
the relative ages docs not lead to an increase in reduced 
(x^), consistent with the scatter in relative ages aris- 
ing from the formal errors alone. Note that freezing the 
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ages introduces no changes to the correlations of A[Z/H], 
A[a/Fe], or A[a/N] with logcr discussed in the previous 
section. This is shown graphically in Figure Oa). Note, 
too, that the mean per degree of freedom has actu- 
ally decreased, providing further proof that fitting for 
the ages does not contribute any statistically significant 
information about the stellar populations. 

Since the relative ages can be constrained without af- 
fecting the other parameters, we now explore whether 
A[a/C] and A[a/Ca] have any significant variation 
among the E/SOs in CL1358-H62. Figures Cfb) and (c) 
show that these additional variables do not improve the 
quality of the fits to the blue Lick indices, and that the 
carbon and calcium abundance ratios show no statisti- 
cally significant variation or correlation with logcr. We 
also verified that the inclusion of these abundance ratios 
has no effect on the correlations of A[Z/H], A[a/Fe], and 
A[a/N] with velocity dispersion shown above. 

Given the extreme sensitivity of the C4668 index to 
both [Z/H] and [a/C], we find these results quite re- 
markable. First, the lack of a correlation of the car- 
bon abundance ratio with velocity dispersion indicates 
that the broadening corrections were accurate (fitting the 
uncorrected line strengths leads to an anti-correlation 
of A [a/C] with logtr). In the context of the steep 
correlation of nitrogen abundance with velocity disper- 
sion, constant carbon abundance ratios have impor- 
tant ramifications, to be discussed in the next section. 
Lastly, though the constraints on the calcium abun- 
dances are not strong, we see no evidence for a trend 
of A[a/Ca] (or, by extension, A[Ca/Fe]) with veloc- 
ity dispersion, contr ary to what has been published 
by p revious authors llSaglia et al.l l2002t iCenarro et al.l 
20031 iF alcon- Barroso et alJl2 003: Michicl sen et alJl2003l: 
Thomas. Maraston. fc Bendeni2003bi|) . 

5. SUMMARY OF THE MODEL FITTING AND 
IMPLICATIONS 

The models explored in this paper indicate: (1) 
that the ages of the elliptical and lenticular galaxies 
in CL1358-I-62 are remarkably homogeneous, consistent 
wit h the analys is of the fundamental plane of these galax- 
ies l|Kelson et al . 2000c); and (2) the stellar populations 
in these galaxies obey tight, well-defined [Z/H]-, and 
[a/N]-tT relations. Together these trends, when propa- 
gated through stellar population models, should repro- 
duce all of the observed early-type galaxy scaling rela- 
tions in CL1358-I-62 and to this end we show the observed 
line strength-line width relations for the sample in Figure 
|S| Note that the figure also includes those galaxies that 
did not satisfy our criterion of having measurements for 
all of the indices fit in the previous discussions, showing 
a total of 33 E/SOs. This full sample is inhomogeneous 
in the sense that any galaxy's stellar population parame- 
ters would not necessarily derive from an identical set of 
indices. Note that along with the eight fitted indices we 
have also included H/3, Fe5015A, and Mg& in the figure. 

5.1. The Line Strength-Line Width Projections 

In the individual line strength-line width diagrams we 
project the metallicity- and nitrogen-a relations, assum- 
ing no correlation of a/Fe with velocity dispersion. The 
lines show loci of constant ages ranging from 1 Gyr 
to 7 Gyr (the latter essentially the defined age for the 



reference point used in recalibrating the models to the 
data; see i i3.1(l . The oldest isochrones clearly reproduce 
the observed line strength-line width relations well in 
CL1358-1-62, and it is worth noting that the slopes of the 
isochrones are predicted to evolve with age. This "evolu- 
tion" arises from the nonlinearity of the model grids. The 
implication is that any observed evolution in the slopes 
of line strength-line width relations with redshift should 
not necessarily be taken as proof that stellar population 
ages vary along the sequence of early-type galaxies. 

Previous analyses of redder Lick indices explicitly re- 
lied on Mg line strengths, and strong correlations of 
[a/Fe] with veloci ty dispersion ha ve been inferred from 
Mg2 or Mgb (e.g. I^rievich et al.lll98lt lGonzalezllT993 
Gorgas et al. 1993; Trager et al. "199^; Kuntschncr et alj 
M)l; Sanchcz-Blazquez 2004; Thomas et al. 2005j). 
However, our data simply do not reproduce the long- 
standing correlation of a-enhancement with a. The 
blue line strength- line width relations of the E/SOs in 
CL13562 are fit well with no correlation, or perhaps 
even with a mild anti-correlation [a/Fe] with logcr. Fur- 
thermore, we note that despite explicitly excluding H/3, 
Fe5015, and Mg& from the fitting, their line strength-line 
width relations are well described by the projections of 
the abundance-velocity dispersion relations found with 
the bluer indices (though our Mgfe observations have 
larger scatter due to uncertainties in the corrections for 
telluric absorption). 

Our estimates of a/Fe have implication s for the origin 
of th e Mg2-, Mg6-, and (Fe)-cr relations ijThomas et alJ 
120051 e.g.). Using the models we can recast the 
metallicity- and nitrogen-cr correlations as Mg-cr rela- 
tions and test whether our results are consistent with 
previous measurements of its slope (though the nitro- 
gen has no impact on the Mg indices in the current 
incarnations of the models). After passively evolving 
these relations to the present epoch, we find Mg2 oc 
(0.23 ± 0.05) logcr, in excellent agreement with that 
found nearby (typi cal slopes nearby range fro m 0.17 
to 0.23 mag/dex; Worthev & Collobert' 2003"). For 
Mg:6, the T homas. Maraston. & Bender (2003a) and 
iTh omas. Maraston. fc KornI lj2004|) models predict that 
Mg6 is about twice as sensitive to changes in a/Fe, so a 
comparison to Mg6 should provide a stronger test of the 
consistency of our results with previously published data. 
We infer Mg6 oc (2.8 ± 0.7) log a , consistent with the 
Mg6 oc 2.8 logcr typically seen fsee lWorthev fc Colloberd 
120031) . Fitting the Mg6 data of the full inhomogeneous 
sample directly, we obtain Mg6 oc (2.1 ± 0.8) logcr (see 
Fig. ^. 

Because our results are consistent with the Mg6-cr (and 
Mg2-cr) relation, we now explore potential sources of the 
discrepancy over a/Fe. In Figure|^our [Z/H]- and [ct/N]- 
cr correlations have been propagated through the models 
and plott ed over data taken in a survey of nearby clus- 
ters from iNelan et al.l l|2005D . Those authors observed 
more than 4000 red-sequence galaxies in 93 low-redshift 
clusters. Because of the large number of galaxies in their 
survey, and the low typical signal-to-noise ratio of their 
spectra, we bin their data in intervals of 0.05 dex, and 
use only those galaxies with minimal detected QUI emis- 
sion. The vertical bars mark the bounds within which 
90% of the galaxies fall. The points show the median 
line strength within each bin. 
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Fig. 8. — The line strength-line width relations for E/SOs in CL1358-I-62. Symbols are as in Figure|S] The solid lines show curves of 
constant age and constant a/Fe, defined by projecting the A[Z/H]-, and A[a/N]-(T relations. The ages range from 1 Gyr to 7 Gyr (violet, 
blue, green, orange, etc.) in increments of 1 Gyr. Several items are worth noting: (1) while H/3, Fe5015, and Mgb were not used in the 
fitting, the fitted correlations of stellar population parameters with velocity dispersion reproduce the observed H/3-, Fe5015-, and Mgb-cr 
relations in CL1358-I-62 (though Mgb has some residual contamination of telluric absorption); (2) we also show all the E/SOs in the sample, 
including those galaxies that did not have observations of all of the fitted indices and were not included in the analysis of the homogeneous 
set of 19 E/SOs discussed in the text; and (3) the slopes of the line strength-line width relations arc sensitive to the ages of t he galaxie s, even 
when there is no variation of age with velocity dispersion. The Thomas, Maraston. & Bender 1 2003ah (and IThomas. Ma raston. &: KornI 
120041) stellar population models do an excellent job of describing the relative properties of the E/SOs in this cluster. 



Overall the comparison between the lNelan et al.l l)2005fl 
data and the predicted line strength-line width relations 
is quite good. However there are particularly interest- 
ing discrepancies: the relations for Fe4383, Fe5015, and 
(Fe), Fe5406, Fe5709, Fe5782. For these indices, our best- 
fit stellar popul ation trends produc e steeper slopes than 
observed in the iNelan et al.l ||200,'tD data. Because (Fe) 
has historically been one of the primary constraints on 
metallicity, we c onclude that the shallow slope of the 
(Fe)-(T relation of INelan et al . (2005) is characteristic of 
the kinds of relations that have led to flatter [Fe/H]-cr 
relations than indicated by our data. Such a flat [Fe/H]- 
(T relation invariably necessitates a positive correlation 
of [a/Fe] with velocity dispersion in order to match the 
slope of the Mg6-cr relation. As can be seen in Fig. |H| 
our Fe4383-(T relation (and Fe5015-a, though with larger 
scatter) has a slope that agrees with the projection of the 
models, thus, in tandem with the other indices, negat- 
ing the need for a significant correlation of a/Fe with 
velocity dispersion. Clearly the primary source of the 
discrepancy in a/Fe comes from the narrow Fe indices. 
Note, too, the divergence of Ca4227 at large values of a. 

These discrepant indices are some of the narrowest 



in the Lick system, and because the remaining indices 
(those that require very small corrections) are matched 
very well by our abundance patterns, one must sus- 
pect that the discrepancy in a/Fe-a arises from differ- 
ences in the treatment of the corrections for instrumen- 
tal and Doppler broadening. These narrow indices have 
had large, systematic co rrections for Doppl er broaden- 
ing applied to them by INelan et alJ l)2005t) (see, e.g., 
Ppggianti et al. 2003, for plots of such corrections as 
functions of velocity dispersion). Systematic errors in 
such corrections can modify the inferred slopes of line 
strength-line width relations. 

We tested whether the different a/Fe-a correlation was 
an inherent property of the sample or whether the treat- 
ment of the corrections played a role. For the first test we 
derived line strengths using the traditional method, i.e. 
first smoothing the spectra to the resolution of the IDS, 
then measuring the indices, and finally applying multi- 
plicative corrections to these measurements. These ve- 
locity dispersio n corrections were constru cted using an 
old, metal-rich iBriizual fc CharlotI l|200l model SED. 
These line strengths were propagated through the ma- 
chinery to derive new stellar population parameters and 
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Fig. 9. — A comparison of isochrones following the [Z/H]- and [a/NJ-cr relations of the text with selected data from the survey of nearby 
clusters by Nelan ct al. (2005;). The points show the median line strengths of those galaxies with OIII equivalent widths less than the 
median of the sample, binned in velocity dispersion in units of 0.05 dex. The vertical bars show the bounds within which 90% of the sample 
are distributed. The predicted loci of galaxies agree well with the distribution of galaxies in the nearby clusters for many of the indices. 
The data of Nelan ot al. (20(3) show mildly flatter line strength-line width relations for the narrow atomic indices Fe4383, Fe5015, (Fe), 
Fe5406, Fe5709, and Fe5782 in contrast to the excellent match to Fe4383 in Fig. |H] These discrepancies, as well as the divergence of Ca4227 
at large dispersion, may be due to differences in the treatment of the broadening corrections as these are the narrowest indices. 



the results are compared in Figure ITUI to the parameters 
found earlier. As can be seen in Fig. Iior b). we recovered 
A[Q;/Fe] oc (-|-0.27±0.13) logcr, in much better agreement 
with previously published correlations. In this test, the 
corrections were computed using a single template. Us- 
ing the best-fit Bruzual & Chariot ( 2003) SED for each 
galaxy instead of a single template, we find no significant 
change in the slope of the inferred correlation (though the 
scatter is larger). 

Given this discrepancy, we performed a second test 
of our broadening corrections by smoothing the galaxy 
spectra to have intrinsic broadenings equivalent to the 
Lick/IDS resolution with zero velocity dispersion, i.e. by 
subtracting, in quadrature, the velocity dispersion from 
the IDS resolution. We remeasured the indices in these 



spectra and compared them to our broadening-corrected 
indices derived earlier. We found no statistically sig- 
nificant systematic differences between these measure- 
ments and the measurements corrected in ti2.4l for the 
line strengths of those galaxies with velocity dispersions 
less than the IDS resolution at the wavelength of a given 
central passband. Therefore we conclude that our treat- 
ment of the broadening corrections has not introduced 
any gross systematic errors as a function of velocity dis- 
persion. 

Lastly, it is worth re-emphasizing that our diagnostics 
of a-enhancements differ from those traditionally used. 
For example, our spectral coverage does not provide di- 
rect constraints on magnesium abundances. The blue 
indices should have stronger (though indirect) sensitivi- 
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Fig. 10. — A comparison of the best-fit parameters Alogi, A[Z/H], A[a/Fe] (wit h no additional nitrogen enhancement) using (a) the line 
strengths corrected for instrumental and Doppler broadenings using the method in i|2.4l (b) the line strengths determined using a traditional 
approach to correcting the effects of instrumental and Doppler broadenings. Subtle differences in the treatment of the broadening corrections 
lead to different inferred correlations of a-enhancement with velocity dispersion. In (a), the best-fit slope is [a/Fe] oc (—0.01 ± 0.18) logcr, 
while in (b), the best-fit slope is [o/Fe] oc (+0.27 ± 0.13) logcr, more consistent with previously published results. 




0.9 



100 200 
a (km/s) 




0.8 



0.9 



(B-V) . 

^ ^ obs 



> 




100 



200 
a (km/s) 



Fig. 11. — A comparison o f observed restframe (B — V)obscrvcd colors iKelson et al.ll2'000al l3) with those predicted using the best-fit 
stellar population parameters IBruzual fc Ch arlot"2003l) . The top panels show predicted colors derived when age was kept constant. The 
bottom panels show the predicted colors when age was allowed to vary. In the left-hand panels the black solid line indicates a slope 
of unity. When only the abundances were allowed to vary, the best-fit slope is (B — ^)j2/H] (0.97 ± 0.11){B — V)obscrvodi shown 

by the green solid line (where the 1-cr uncertainties are shown using the green dashed lines). When age is allowed to vary one obtains 
{B -— V)^ [Z/H] °'~ (^'^"^ -'- 0-16)(i3 — V)obscrvodi also shown in green. The effects of including age are small because of (1) the statistically 

insignificant variations in age in the sample, and (2) the large covariances between [Z/H] and logt and the age-metallicity degeneracy in 
{B — V) colors. The differences between the observed and predicted colors are shown plotted against velocity dispersion in the right-hand 
panels. The residuals in the colors are remarkably small, consistent with the formal uncertainties. 



ties to oxygen abundances, primarily through the CNO 
equilibrium in the atmospheres of cool giants. Because 
oxygen is the dominant a element, our absence of a cor- 
relation of a/Fe with velocity dispersion may be a more 
accurate descriptor of the behavior of a elements than 
indicated by Mg alone (though we note the excellent 



match to the slopes of published Mg-cr relations). In 
other words, our a/Fe results may simply reflect more 
on a uniformity in 0/Fe than in Mg/Fe, not to mention 
our empirically constrained uniformity in the carbon and 
calcium abundance ratios. 
While the broadening corrections appear to play an 
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important role in the determinations of a/Fe in the 3- 
parameter fitting, these corrections become less impor- 
tant when including the nitrogen enhancement. When 
using the line strengths corrected for broadening via the 
traditional method, we obtain similar results to those 
shown in Figure Et^c) for the 4-parameter fitting of age, 
metallicity, a/Fe, and N/a. In ti5.3l the inferred nitrogen- 
a correlation will be discussed in more detail, coming to 
the conclusion that it has a natural origin. 

In the next section we show that the observed colors 
of the galaxies match those predicted by the stellar pop- 
ulation parameters that we have derived. 

5.2. The Color-Magnitude Relation 

The [Z/H] -cr correlation has strong implications for the 
color-magnitude relation. We e xplore these using the 
WFP C2 data pu blished by |^ Dokkum ct al. ( 1998a) 
and Kelso n et all (pOOOa). These images were convolved 
with Gaussians to mimic the seeing in our spectroscopy. 
Then (F606W— F814W) colors were measured within cir- 
cular apertures similar to the apertures used for the spec- 
tral extractions. The colors were transformed to rest- 
frame Johnson colors, {B — V^) observed- In this section we 
compare these observed c olors to those predicted by the 
iBruzual fc Chariot l)2003|) stellar population models us- 
ing ages and metallicities from our analysis. In order to 
eliminate any systematic errors in the photometric trans- 
formation, in the calibration, or in the uncertainties in 
the absolute ages and metallicities, we applied an offset 
to the colors computed us i ng the biweight estimator of 
IBeers. Flvnn. &: Gebhard3 ljl99Cf) . much like what was 
done earlier for the line strength fitting. 

FigurelTTlshows a comparison of these predicted colors 
with the observed {B — y)observed colors. When the ages 
were kept constant, we refer to the predicted colors as 
{B — V) \^ ^jjj . When the ages were allowed to vary in the 

fit, the predicted colors are referred to as {B — V^)^ [Z/H] • 
In the left-hand panels one finds that both sets of pre- 
dicted colors match the observed colors well. Simple fits 
to the comparisons yield (S— l^)[z/H] (0.97±0.11)(i3— 

^) observed ; 

and {B - ^)t,[z/H] « (1-14 ± 0.16)(i? - 
V')observed, and thcsc are shown using the green lines. 
The solid black line indicates a slope of unity. Because 
colors suffer from an age-metallicity degeneracy similar 
to that which plagues the absorption line strengths, the 
comparison of the colors only serves as a sanity check 
— stellar populations are indeed redder as they become 
older or more metal-rich, and in the correct proportion 
to the changes in spectroscopic ages or metallicities. In 
other words, one gains no additional leverage on the stel- 
lar population parameters by fitting (B — V) colors with 
the absorption line strengths. 

Simply put, the fitted ages and metallicities repro- 
duce the color-velocity dispersion relations remarkably 
well. In the right-hand panels we plot the color resid- 
uals against a, and the very low scatter, ±0.02 mag, is 
consistent with the formal uncertainties. 

Had the ages and metallicities from Figure [TUT b') been 
used, the slopes of the correlations between the pre- 
dicted and observed colors would have resulted in slopes 
of (S - l^)[z/H] « (1-25 ± 0.10)(B - y)observed, and 

{B - ^)t,[Z/H] « (1-50 ± 0.17)(B - y)observed, with sig- 



nificantly larger scatter. This is because a significant, 
positive correlation of a-enhancement with velocity dis- 
persion decreases the slope of the [Z/HJ-cr correlation to 
the point that the color-velocity dispersion relation can 
no longer be matched. 

In summary, the colors predicted by our population 
parameters are in excellent agreement with the observed 
colors, when measured in the appropriate aperture. This 
remarkable consistency between the colors and the spec- 
troscopic properties of these galaxies also implies that 
dust is not an important contributor to the colors of 
early- type galaxies. 

5.3. The Nitrogen Relation 

In fitting the iThomas. Maraston. fc Benderf l)2003aj) 
models to the blue line strengths, we found a strong ap- 
parent correlation of nitrogen enhancement with velocity 
dispersion. However, there appears to be no associated 
correlation (or anti-correlation) with the carbon abun- 
dance ratio. The lack of concordance between the nitro- 
gen and carbon enhancements was not expected. Taken 
together with our inability to reproduce the previously 
published correlation between [a/Fe] and log cr, our re- 
sults would seem to indicate either fundamental prob- 
lems with the stellar population models, or something 
truly fundamental about the formation and enrichment 
histories of galaxies. 

Certain ly, uncertainties in the models are important 
(see e.g. iGratton. Sneden. &: Carettal 12004 for a re- 
view of unexplained C and N variations among glob- 
ular clusters). The nitrogen- and carbon-sensitive 
indices, are affected by mixing from first dredge- 
up, by other mixing processes as well (again, see 
IGratton. Sneden. fc Carettal 12004 for a thorough sum- 
mary and detailed references), and by the unknown 
abundances of oxygen, which are a wholly s eparate 
canard on their own (Bcnsbv. Feltzing. & Lun dstron] 
12004 IProctor et al.ll2064: .Thomas et al...2005,) . Recent 
work on understanding the carbon and nitrogen abun- 
dances in Galactic and globular clusters continue to 
confound (e. g. Aoki & Tsuii 1997; Bellman ct al. 200Ji 
iBrilev. Cohem fc StetsonI I2002L l26o4a(l . though it ap- 
pears likely that both pollution by early generations of 
intermediate-mass a symptotic giant branch (AGB) stars 
(jBrilev et a l.' '2004b'. e.g.) and additional mixing pro- 
cesses (e.g. .Boothrovd & Sack mann 1999) occur. Taken 
all together, one cannot rule out if the deduced N/C 
correlation with velocity dispersion is an artifact aris- 
ing from an absence of processes that modify the atmo- 
spheric abundance ratios in the models, processes that 
arc themselves sen sitive to metalli city (e.g. Charbonne| 
^994; Bens bv. Feltzing. fc Lundstromil2004|) . 

Recalibration of the models to the data eliminated 
several systematic errors, as illustrated in Appendix 
1X1 and as a result, the relative measures of age 
and abundance are probably quite reasonable. How- 
ever, the N/C ratios may yet be reflecting limi- 
t at ions in the s tellar libraries and fitting functions 
llWorthevI 179941 IThomas. MaTa,ston. fc BendeH 12003^ 
iThomas. Maraston. fc Kornll2004^ . in that the param- 
eter space of second derivatives of the indices has not yet 
been fu lly expl ored and modeled iTripicco fc BeU 19951 
iKorn. Maraston. fc Thomas! 120051) . Sufficiently large, 
but unmodeled, non-linearities in the underlying sensi- 
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Fig. 12. — The correlation of A[a/N] and A[Z/H] . For simple models of chemical enrichment IKoppen fc Edmundsl 
IHenrv. Edmunds. Hi KopDenll2000l : IKoppen fc Hensleiti2005l) in which abundances are dominated by secondary nitrogen, one expects the 
relation to follow A[a/N] = — A[Z/H], shown by the solid line. We fit for the best-fit linear relation, using those galaxies with a > 134 
km/s, by minimizing the squares of the residuals parallel to both axes of the error ellipses simultaneously. The resulting best-fit slope is 
(— 1.07±0.10) dox/dex, shown by the dashed and dotted lines. The agreement between the stellar population parameters and the prediction 
of simple chemical evolutionary models is striking. Note that this slope was derived using galaxies above the limit of the selection, and 
thus excluded, for example, galaxy 493 with the lowest metallicity and highest nitrogen abundance. 



tivities of the indices could be invoked to eliminate the 
nitrogen abundance trends but the efforts to calibrate 
the mo dels using the integrated light of globular clusters 
(e.g. Schiavon et al. 2002: T homas. Maraston. fc Bendeil 
l2C)03a: Proctor et al. 2004) provide some confidence in 
our results. 

It is worth noting that the more common line strength- 
line width relations, such as Mg-cr, appear to be well- 
described without any need for a correlation of a- 
enhancement with velocity dispersion, all while the blue 
Lick indices are fit quite well. This self-consistency did, 
however, come at a cost: that of throwing away the zero- 
points of the models and the absolute measures of abun- 
dance and age. Had we not done so, the reduced on 
individual galaxies increase by factors of 2-7! Taken to- 
gether with the large systematic uncertainties shown in 
the appendix, we have much greater confidence in the 
relative measures of abundance derived above than in 
estimates of absolute abundances. 

The steep correlation of the nitrogen abundance ra- 
tio with velocity dispersion has interesting implica- 
tions, but perhaps not for the most obvious reason. 
Naively one might invoke a physical correlation between 
gravitational potential and the timescales and/or effi- 
ciency of star-form ation (e.g. Hcnrv & Worthcv 1999; 
Im)rthevl ll998: He nrv. Edmunds, fc Kopp en 2000) given 
the sensitivities of nitrogen abundances to these pa- 
rameters for closed-box models. However, the under- 
lying correlation for nitrogen is probably not with ve- 
locity dispersion per se but with mean metallicity: sec- 
ondary nit rogen is produced in intermediate-m ass AGB 
stars (e.g. IHenrv. Edmunds, fc Koppenl 120001) and th e 
yields are metallicitv-dependent ijMarieo et aP Il998j) . 
Closed-box (and partially closed-box) models of chem- 
ical enrichment with modest or high metallicities have 
N/0 ratios dominated by s econdary nitrogen such that 
A[N/01 (X A[0/H] (e.g. IKoppen fc Edmundsl IT99i 



Henrv. Edmunds, fc Koppen 1200(1 iKoppen fc Henslen 
2005ft . The dependencv of the yields on metallicity sets 
the absolute zero-point of the logarithmic relation, but 
one expects the relation between relative nitrogen abun- 
dance ratios and mean metallicities to have unity slope. 



From the tables in IHenrv. Edmunds, fc K5ppe n (2000), 
the carbon yields of AGB stars do not appear to be par- 
ticularly sensitive to metallicity; presumably this is why 
our carbon abundance ratios do not correlate with ve- 
locity dispersion (or metallicity) in the same way that 
nitrogen does. 

The hypothesis that the a/N-cr relation is due to sec- 
ondary nitrogen can be tested because [N/0] is predicted 
to correlate with [O /H] . Within the context of our data, 
[a/N] is simply a surrogate for [0/N]. Furthermore, even 
though oxygen abundances cannot be directly measured, 
our findings of A[a/Fe] ~ (or A[0/Fe] ~ 0) imply that 
we can approximate A[Z/H] ~ A[0/H]. In Figure IT^ 
the relative nitrogen abundance ratios A[a/N] against 
A[Z/H]. Finding the best-fit relation between A[a/N] 
and A[Z/H] is complicated by the correlated errors in the 
derived abundances so we chose to minimize the squares 
of the residuals parallel to both axes of the error ellipses 
simultaneously. We find that the slope of the correlation 
is (—1.07 ±0.10) dex/dex, shown by the dashed and dot- 
ted lines in the figure. Using other approaches to fit the 
correlation led to similar results. We therefore conclude 
that our data are consistent with the hypothesis that the 
[a/N]-CT correlation arises from the dependence of yields 
of secondary nitrogen on metallicity. 

The stellar populations of the E/SO galaxies in 
CL1358-I-62 have essentially been reduced to a one- 
parameter family, obeying a [Z/H]-(t relation. While 
the presence of (secondary) nitrogen in proportion to the 
metallicity indicates that the galaxies experienced star- 
formation over an extended period of time, absolute cali- 
bration of the abundances are required in order to derive 
constraints on the timescales of astration. Even without 
absolute calibrations, it bears repeating: the populations 
of both Icnticulars and ellipticals have statistically identi- 
cal relative ages, span similar ranges of metallicity, obey 
identical [Z/H]-(t relations, and show strong evidence for 
the presence of secondary nitrogen. 

Our results conclusively show that cluster SOs did not 
form their stars at significantly later epochs than clus- 
ter ellipticals of the same mass, and the presence of sec- 
ondary nitrogen indicates that both Es and SOs formed 
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from ancestors with extended star-formation histories. 
If our a/Fe results are correct, then this star- formation 
must have been sufficiently extended, within a given 
galaxy, such that the relative proportions of by-products 
of Type la to Type 11 supernovae are effectively con- 
stant, and not strongly dependent on galaxy mass. Fur- 
thermore, if the gravitational potential is an important 
factor in the retention or ejection of enriched materi- 
als (e.g. I&arke 19(3 and much subsequent work), then 
material from all energetic sources of enrichment have 
roughly comparable efficiencies of ejection from the host 
galaxies, either because of generally low probabilities of 
ejection, or because the enriched interstellar gas is well- 
mixed before it either forms the next generations of stars 
or escapes from its host. 

This produces an apparent contradiction with the 
obser ved evolution i n th e morphology-density rela- 
tion ijPressler et al.l Il997t iPostman et ail 12005(1 . in 
which fewer SOs e xist at higher redshift (though see 
IHoldeu et a,].l [200fih. even though the timescales and 
mechanisms with which SOs form must be very similar 
to the timescales mechanisms for the formation of ellip- 
ticals of the same mass. 

Ultimately, the chemical signatures, and any evolution 
thereof, must be placed in the context of cosmological 
models. In ACDM, galaxies form hierarchically, th rough 
multiple generations of r nergers (. Kauffmann et al JIlQQJ^ 
ir999- .Romeo et al.ll2005l iLanzoni et al..i2005t .Tran et all 
120051) . Unfortunately, multiple major merger events 
can have a lar ge impact on the observed scaling rela- 
tions fe.g. lKavir ai et al. 2005), such as the line strength- 
line width relation s analyzed here. Furthermore, (e.g. 
I.Tuneaii et al.l2005l) directly measured the star-formation 
rates of galaxies with these masses over redshifts of 
0.8 < z < 2 and found that their star-formation has 
essentially ceased by z > 1.5. In a subsequent paper we 
will try to address these issues by analyzing the stellar 
populations of both early- and late-type galaxies in a se- 
quence of five clusters from z — to z — 0.83, referencing 
them to the galaxies presented in this paper. Together it 
is hoped that these five clusters can provide a more de- 
tailed picture of the formation histories of cluster galaxies 
than can be derived from a snapshot at a single epoch. 

6. CONCLUSIONS 

We have measured absorption line strengths for galax- 
ies in th e cluster CL 1 358-I-62 using the spectra pub- 
lished in iKelson et"an l(2000bl ^. A homo geneous popu- 
lation of early-type galaxies has been ana l yzed, w ith the 
full selection of galaxies in IKelson et al.l l|2000b() to be 
discussed in a subsequent paper. The largest source of 
systematic errors that plague the analysis of absorption 
lines has been eli minated b y matchi ng the stellar pop- 
ulations models ||Thonias^Marasto n. fc Bendeil l2003al: 
iThomas. Maraston. fc Kornll20041 to the mean observed 
line strengths of the most massive early-type galaxies 
in CL1358-I-62. Furthermore, not only does a recalibra- 
tion of the models allow for accurate differential mea- 
surements of the stellar population parameters, but the 
topology of the minimum can be used to accurately 
estimate formal errors. 

Using only those E/SO galaxies which have accurate 
measurements of eight blue Lick/IDS indices (11(5^, Il7^, 
CNa, Ca4227, G4300, Fe4383, Fe4531, and C4668), we 



fit for the relative measures of age and chemical abun- 
dances. The fitting has resulted in homogeneous sets 
of stellar population parameters relative to the reference 
point used in resetting the model zero-point. 
The key results can be summarized as follows: 

(1) The populations of the E and SO galaxies are sta- 
tistically identical in their age and abundance patterns, 
down to the magnitude limit of i? = 21 mag, consis- 
tent with what was found using the fundamental plane 
(|Kelson et al.ll20'00d) . Beyond the magnitude limit, se- 
lection biases become important. 

(2) The observed scatter in relative ages is 0.06 
dex, where the scatter expected from measurement er- 
rors is 0.09 dex. The scatter in ages is consistent 
with t hat inferred from the color-magnitude relation 
l(van D okkum et al. 1998as), and from the fundamental 
plane (Kelson et al. 20003). 

(3) We find a tight correlation between relative metallic- 
ity and velocity dispersion, with a slope of 0.86 ± 0.17 
dex/dex. The scatter about this relation is 0.06 dex, 
where the scatter expected from the formal errors is 0.05 
dex. The observed B — V colors are also consistent with 
the inferred ages and metallicities. 

(4) The scatter in relative a-enhancement is only 30% 
larger than expected from the formal errors alone. When 
the nitrogen abundance ratio is allowed to vary, we find 
a mild anti-correlation of A [a/Fe] with logcr, but with 
a significance less than < 2a. While we find little evi- 
dence for a significant correlation (or variation) of [a/Fe] 
with velocity dispersion, the slope of the local Mg-cr re- 
lation ca n be reproduced by our A[Z/H]-(t correlation. 
iThomas et al.l (2005) argued that Mgfe-cr arises mostly 
from the correlation of metallicity with velocity disper- 
sion. Our results go beyond that and suggest that no 
variation of a/Fe is required. The discrepancy between 
our conclusions and previously published work on this 
topic arises from our steeper [Z/Hj-cr relation. This dis- 
crepancy may be due to the treatment of the broadening 
corrections, in particular for the narrow Lick iron indices, 
though more detailed analysis is required to verify this 
hypothesis. Previous work (e.g. Trager et al. 2000b) in- 
dicated that cluster E/SOs form a three-parameter family 
of objects, with [Z/H] and a/Fe correlated with velocity 
dispersion. Historically, models with variable ratios of 
Type I and Type II supernovae and subsequent galac- 
tic winds have been invoked to explain such a complex 
family of objects. Elimination of the correlation of a- 
enhancement with velocity dispersion would greatly sim- 
plify models of E/SO formation. 

(5) We find that the enhancement of nitrogen is strongly 
correlated with velocity dispersion, while carbon abun- 
dance appears to show no variation. This N-ci correlation 
has not been observed before. This new correlation origi- 
nates either in deficiencies in the stellar population mod- 
els, or implies new physics to be incorporated into our un- 
derstanding of E/SO formation and evolution. The most 
likely explanation is one in which the nitrogen abundance 
ratios do not arise from mechanisms specifically tied to 
tr, but that the nitrogen-enhancement is correlated with 
metallicity. In other words, the nitrogen is secondary in 
origin, and the progenitors of E/SO galaxies experienced 
significantly extended star-formation. 

These data form the only sample of early-type galaxies 
in which these particular abundance patterns have been 
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found, so far. In all other respects the sample shows the 
uniform age distribution and a steep correlation of mean 
metallicity with velocity dispersion inferred from other 
diagnostics. 

W hile models w hich inclu de variable abu ndance ratios 
(e.g. [ T rager et aDr2 000b: T homas. Marast on. fc Bender 
I2003al iThomas. Ma raston. fc Kornll200^ are relatively 
immature, their development has been encouraging, and 
the data are already well-fit by such models. The blue 
Lick indices appear to be powerful diagnostics of stel- 
lar population ages and provide excellent, self-consistent 
fits to the metal lines, even if the resulting abundance 
ratios run counter to much of the previous literature. 
Despite these discrepancies, the observed line strength- 
line w idth relations of nearby galaxies, published by 
INelan e t al. (2005), are recovered, with some notable ex- 
ceptions. Those relations that match well do bolster our 
faith in these findings. However, nearly all of the iron 
line strength-line width relations are not matched well, 
and the published data on nearby galaxies have had large 
multiplicative Doppler corrections applied. We suggest 
that perhaps past treatment of these corrections may be 
in error. That the {B — V) colors are consistent with our 
stellar population parameters not only reinforces our con- 
clusions, but also indicates that dust is not an important 
component in the optical SEDs of early- type galaxies. 

The results on enhancements (over Fe) of a, nitro- 
gen, and carbon have important implications. While it 
remains to be seen whether the carbon- and nitrogen- 
sensitive indices are being modeled correctly, the nitro- 
gen enhancement-velocity dispersion correlation agrees 
very well with the correlation of [N/Z] with [Z/H] ex- 
pected to arise from secondary production of nitrogen 
by AGB stars. Taken together, our differential stellar 
population parameters suggest that these massive early- 
type galaxies form a family with only a single parameter: 
galaxy velocity dispersion (or possible mass). 



The uniform ages and simple trends of metallicity with 
velocity dispersion appear to contradict t he observation 
that f ewer SOs exis t at higher redshift ijDressler et alJ 
Il997t iPostman et~ l. 2005). We conclude that the 
timescales and mechanisms with which the populations 
in SOs form are statistically identical to that of ellipti- 
cals of the same mass. Early- type galaxies, now shown 
to have significant self-enrichment and extended star- 
formation histories , must be reconciled w ith hierarchi- 
cal formation (e.g. Kauffm ann et all 11993!) in a ACDM 
universe. The uniform ages and tight mass-metallicity 
relation may not be compatible with significant merging 
and histories of extended star-formation. 

In a subsequent paper, we will employ our machin- 
ery for measuring these age and abundance patterns in 
a range of clusters from z = to z = 0.83. With 
comparisons of the properties of cluster galaxies over a 
range of redshifts, we will test the universality of some of 
the conclusions presented here, quantify changes in the 
stellar population parameters with cluster redshif t (see, 
e.g., K elson ct al. 2001; J0rgcnsen et al. 2005; Bar r et alJ 
|2005|) . and attempt to reconcile the above results with a 
coherent picture of the formation and evolution of galax- 
ies in rich clusters. 
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APPENDIX 

THE EFFECTS OF THE SYSTEMATIC RECALIBRATION 

Differential versus Absolute Population Parameters 

In this section we show several figures to illustrate the effects of recalibrating the models to our data. In particular 
we comp are fig ures shown earlier to variants in which this step is ignored. 

Figure IA13I shows the parameters of t he be st-fit models, in which age, metallicit y, g- enhancement, and nitrogen 
enhancement are allowed to vary. Figure lA13r a) is identical to Fig. IHI^c). In Figure lAlSf b) we show the parameters 
derived when the zero-point uncertainties are ignored and the models and data are assumed to be on the same system. 
Similar trends are seen in both the differential and absolute parameters. However the absolute parameters, derived 
when the recalibration is skipped, have 30-40% larger scatter. These systematic uncertainties appear not to destroy 
the general trends of the SSP parameters with velocity dispersion but artificially infiate the inhomogeneities in the 
inferred properties of the stellar populations. 

In Figure IA14I we show another comparison of the differential and absolute analyses, in particular investigating the 
effects on parameter-parameter correlations. For this discussion we exclude those galaxies below cr = 134 km/s to 
reduce the impact of the luminosity-selection bias. In the left-hand panels we show the relative ages and nitrogen 
enhancements plotted against relative metallicity. No significant correlation exists between the relative ages and 
relative metallicity. The significant correlation between relative nitrogen enhancements and metallicity was discussed 
earlier. If one ignores the recalibration of the models to our data, one obtains an anti-correlation between age and 
metallicity with a slope of —0.76 ± 0.10 de x/dex (shown by the dashed and dotted lines). Note that this artificially 
induced slope obeys the famed "3/4" rule of lWorthevI l| 19941) because the systematic errors in our line strengths couple 
with the average age-metallicity degeneracy of the indices. 

Furthermore, in the absolute analysis the correlation between nitrogen enhancement and metallicity effectively 
disappears. Because Fig. IA13f b) shows the strong [Z/Hj-cr and [N/aJ-cr correlations, a lack of an observed correlation 
between [Z/H] and [N/a] would therefore have masked the true nature of the [N/aJ-cr correlation. Note that while 
the general trends of the SSP parameters with velocity dispersion were preserved when the recalibration was ignored, 
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Fig. A13. — The effects of ignoring the model recalibration on the stellar population parameters as functions of velocity dispersion for 
ellipticals (filled red circles) and E/SO and SO galaxies (open blue circles). Only those galaxies with all eight observables are shown, (a) 
A log t-A[Z/H]-A[a/Fe]-A[a/N] models, (b) log t-[Z/H]-[a/Fe]-[a/N] models. The left-most panels show the reduced for each galaxy, 
also plotted against galaxy velocity dispersion. The mean reduced values are also given in these panels. While the general trends shown 
in Figure|n]are clearly visible, the scatter about the correlations is 30-40% larger when the models are not recalibrated to the data. 



the systematic errors in the cahbration had resulted in larger errors in the SSP parameters. These errors appeared 
in Figure fAlSf b) as larger scatter and increased inhomogeneity. However, these added errors can be both correlated 
and uncorrelated, possibly imposing apparent (anti)correlations between age and metallicity, or erasing underlying 
correlations between metallicity and abundance ratios. 

For samples of galaxies at high redshift, it remains technical challenge to ensure that absorption line strengths from 
the models and data have no systematic offsets. Extreme care must be taken to eliminate all sources of systematic 
error in order to utilize the models' absolute scales of age and abundance. Fortunately, these large systematic sources 
of error can be eliminated if one is only interested in the relative properties of galaxies within an homogeneous sample. 

Uncertainties in the Adopted Model Zero-point 

In ^13.11 we adopted a set of parameters describing the mean stellar population of the most massive E/SO galaxies 
in CL1358-I-62. These parameters are reprinted as Model in Tab le IX5l Here we investigate some of the effects 
one obtains by varying the parameters of this reference point. Table IXHl lists the standard parameters and the three 
additional variations, Models A, B, and C. These variations are not meant as formal uncertainties in age, [Z/H], 
or a-enhancement, but were specifically chosen to be significantly larger than the uncertainties in those parameters. 
As a result, the general effects on our results may be considered very conservative systematic uncertainties in our 
conclusions. Also note that models in which the reference point has solar metallicity and solar a/Fe have been ruled 
out by many a uthors. 

Figure rA_15l shows the best-fit relative stellar population parameters plotted against velocity dispersion, just as in 
Figure Etc). No change in slope in any of the correlations of SSP parameters with velocity dispersion is statistically 
significant. Because the true uncertainties in the zero-point are likely to be much smaller than the range probed by 
our variations in age, metallicity, and a-enhancement, we conclude that our results are not sensitive to systematic 
uncertainties in the recalibration. 
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Fig. A14. — The relative and absolute age-metallicity and nitrogen-metallicity diagrams for those E/SO galaxies with a > 134 km/s. (a) 
and (c) show the results from the differential fitting, (b) and (d) show the parameters derived assuming no zero-point correction. When 
the models are not recalibrated to the data, the results are quite striking: (1) a "statistically" significant anti-correlation between logt and 
[Z/H] appears, with a slope of —0.74 ± 0.13 dex/dex (shown by the dashed and dotted lines); and (2) the nitrogen-metallicity correlation 
completely disappears, masking the true nature of the nitrogen enhancements. Systematic uncertainties in the calibrations of our indices 
would have obscured the true nature of early-type stella r populations. Note that the induced correlation between the absolute ages and 
metallicities obeys the famed "3/4" rule IWorthevlll994ll . 
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Fig. A15. — The effects of varying the model zero-point on the relative stellar population parameters as functions of velocity dispe rsion. 
In (a), (b), (c), and (d) we show the results using the zero-points defined by Models 0, A, B, and C, respectively, defined in Table lASl Note 
that (a) is identical to Fig. IHIc). The general trends, and scatter about these trends, are insensitive to these large changes in the adopted 
zero-point. Because the variations in the recalibration probed here are much larger than the true uncertainties in the reference parameters, 
we conclude that our results are insensitive to systematic errors in the recalibration. 
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Line Strengths in E/SO Galaxies at z=0.33 
TABLE Al 

Corrected Velocity Dispersions and Balmer Line Strengths 



ID 


Type 




log (T 


H5a 




H7A 






H/3 


095 


SO 


2 


321 ± 0.013 


-1.37 ±0.57 


0.22 ± 0.17 


-5.11 ±0.30 


-1.50 ±0.17 


1 


72 ±0.13 


135 


SO 


2 


180 ±0.014 


-0.35 ± 0.28 


1.01 ± 0.18 


-3.83 ± 0.27 


-0.62 ±0.16 


2 


43 ±0.16 


182 


SO 


2 


091 ± 0.018 


-0.35 ± 0.40 


0.38 ± 0.27 


-3.49 ±0.40 


-0.63 ±0.23 


1 


87 ±0.27 


211 


SO 


2 


257 ±0.015 


-0.97 ±0.29 


0.51 ± 0.19 


-4.47 ± 0.31 


-0.93 ±0.19 


2 


21 ± 0.19 


212 


E 


2 


235 ± 0.014 


-0.52 ± 0.28 


0.52 ± 0.20 


-4.41 ± 0.23 


-0.70 ±0.13 


1 


79 ±0.13 


215 


SO 


2 


198 ± 0.015 


-0.88 ± 0.38 


0.50 ± 0.25 


-4.41 ±0..39 


-1.09 ±0.24 


2 


29 ±0.21 


233 


E/SO 


2 


348 ± 0.008 


-0.97 ±0.17 


0.92 ±0.10 


-4.60 ±0.17 


-0.88 ±0.10 


2 


15 ±0.09 


236 


SO 


2 


271 ± 0.012 






-5.10 ±0.37 


-0.83 ±0.18 


1 


86 ± 0.21 


242 


E 


2 


304 ± 0.011 


-1.27 ±0.22 


0.51 ±0.14 


-4.97 ±0.25 


-1.44 ±0.16 


1 


37 ±0.14 


256 


E 


2 


395 ± 0.007 




0.43 ± 0.10 


-4.85 ±0.16 


-1.14 ±0.08 


1 


75 ± 0.08 


269 


E/SO 


2 


482 ± 0.008 


-1.78 ±0.16 


0.21 ±0.10 


-5.45 ±0.20 


-1.37 ±0.11 


1 


84 ± 0.08 


292 


SO 


2 


028 ± 0.012 






-2.21 ±0.25 


0.66 ±0.16 


2 


38 ±0.16 


298 


SO 


2 


449 ± 0.008 


-1.86 ±0.19 


0.51 ±0.12 


-5.29 ±0.21 


-1.39 ±0.13 


1 


85 ± 0.09 


300 


SO 


2 


372 ± 0.010 


-1.29 ± 0.21 


0.61 ±0.15 


-5.21 ±0.26 


-1..50±0.15 


1 


58 ±0.13 


303 


E 


2 


191 ± 0.011 


-0.55 ± 0.22 


0.80 ±0.13 


-4.21 ±0.26 


-0.70 ± 0.15 


1 


94 ±0.15 


309 


E/SO 


2 


336 ± 0.009 


-1.31 ±0.20 


0.37 ±0.13 


-5.49 ±0.24 


-1.51 ± 0.14 


1 


69 ±0.10 


353 


E/SO 


2 


336 ± 0.008 






-5.18 ±0.21 


-1.41 ±0.11 


1 


52 ± 0.09 


359 


SO 


2 


279 ± 0.010 






-4.09 ± 0.25 


-0.73 ±0.15 


1 


92 ± 0.12 


360 


E 


2 


226 ± 0.048 






-3.32 ±0.97 


-1.08 ±0.63 


1 


52 ± 0.52 


375 


E 


2 


470 ± 0.010 






-5.69 ±0.35 


-1.92 ±0.30 


1 


54 ±0.58 


381 


E/SO 


2 


304 ±0.011 






-4.54 ±0.26 


-0.95 ±0.15 


1 


83 ±0.15 


391 


E/SO 


2 


392 ± 0.011 




0.74 ±0.23 


-4.49 ±0.21 


-0.96 ±0.13 


1 


83 ±0.11 


408 


SO 


2 


352 ± 0.011 




0.80 ±0.18 


-4.77 ±0.26 


-1.10 ±0.14 


1 


75 ±0.11 


409 


E 


2 


008 ± 0.027 


-0.85 ±0.51 


0.69 ± 0.29 


-4.39 ±0.58 


-0.62 ±0.34 


2 


39 ± 0.41 


410 


SO 


2 


149 ± 0.015 






-4.71 ±0.41 


-0.99 ±0.23 


2 


54 ± 0.23 


412 


E 


2 


206 ± 0.016 


-0.30 ± 0.27 


0.76 ± 0.18 


-3.85 ±0.27 


-0.56 ±0.18 


1 


52 ± 0.20 


463 


SO 


2 


431 ± 0.010 






-4.52 ±0.22 


-0.97 ±0.13 


1 


51 ±0.14 


481 


SO 


2 


007 ± 0.022 


-0.04 ± 0.42 


1.73 ±0.25 


-2.03 ± 0.48 


0.14 ±0.32 


2 


47 ± 0.31 


493 


E/SO 


2 


043 ± 0.031 


1.85 ±0.57 


1.99 ± 0.38 


-4.38 ± 0.60 


-0.94 ± 0.36 


2 


03 ± 0.44 


531 


E 


2 


426 ± 0.008 


-1.68 ±0.12 


0.35 ± 0.09 


-5.30 ±0.17 


-1.37 ±0.09 


1 


77 ± 0.08 


534 


E 


2 


060 ± 0.022 






-3.51 ±0.43 


-0.35 ± 0.26 


2 


00 ± 0.26 


536 


E 


2 


382 ± 0.009 


-1.68 ±0.20 


0.56 ± 0.13 


-5.28 ±0.22 


-1.45 ±0.13 


1 


79 ±0.10 
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TABLE A2 

Corrected Molecular and Metal Line Strengths 



ID 


CNl 


CN2 


Ca4227A 


G4300A 


Fc4383A 


Ca4455A 


Fc453lA 


C4668A 


Fe5015A 


Ms6 


095 


0.06 ±0.01 


0.08 ± 0.01 


0.69 ± 0.34 


5.11 ±0, 


.26 


4, 


.47 ±0.32 


1 


.04 ±0.19 


3 


.24 ±0.25 


6.02 ± 0.35 


4.83 ± 0.34 


4.27 ±0.29 


135 


0.01 ±0.01 


0.04 ±0.01 


0.75 ±0.15 


5.02 ±0, 


,23 


3 


.73 ±0.39 


1 


.00 ± 0.22 


2 


.97 ±0.28 


3.95 ± 0.76 


5.32 ±0.38 




182 


-0.03 ±0.01 


0.01 ± 0.01 


1.09 ±0.17 


5.12 ±0, 


,36 


2, 


,87 ±0.60 


0, 


,43 ± 0.34 


2 


.81 ± 0.36 


5.07 ±0.57 


2.88 ±0.39 




211 


0.03 ±0.01 


0.06 ± 0.01 


1.27 ± 0.12 


5.35 ±0, 


,28 


4, 


,05 ± 0.44 


1, 


,12 ± 0.36 


2, 


,57 ±0.27 


5.29 ± 0.55 


4.48 ±0.41 




212 


0.03 ±0.01 


0.06 ± 0.01 


0.91 ±0.13 


5.41 ±0, 


,21 


2, 


,98 ±0.34 


0, 


,87 ±0.20 


2, 


,65 ± 0.22 


5.45 ± 0.53 


4.29 ±0.29 


2.74 ±0.38 


215 


0.02 ±0.01 


0.04 ±0.01 


0.78 ±0.19 


4.99 ±0, 


,33 


3, 


,30 ± 0.48 


0, 


,89 ± 0.27 


3, 


,07 ±0.37 


5.68 ± 0.47 


4.38 ±0.36 


2.99 ± 0.43 


233 


0.03 ±0.01 


0.06 ± 0.01 


0.87 ±0.08 


5.61 ±0, 


,13 


3, 


.88 ±0.21 


1, 


.07 ±0.11 


3, 


,34 ±0.16 


5.06 ± 0.26 


4.58 ±0.22 




236 


0.04 ± 0.01 


0.07 ±0.01 


0.96 ±0.17 


5.43 ± 0, 


,26 


4, 


.27 ±0.55 


1, 


.05 ±0.21 


3, 


.39 ± 0.25 


6.14 ±0.51 


4.64 ±0.31 


4.20 ±0.30 


242 


0.05 ± 0.01 


0.08 ±0.01 


1.06 ±0.11 


5.05 ± 0, 


,20 


4, 


.01 ±0.23 


1, 


.18 ±0.15 


2, 


.87 ±0.17 


5.79 ± 0.77 


4.48 ± 0.28 


3.99 ± 0.26 


256 


0.06 ±0.01 


0.09 ± 0.01 


1.08 ±0.11 


5.35 ±0, 


,14 


4, 


.24 ± 0.36 


1, 


.19 ±0.29 


3, 


.31 ±0.14 


6.11 ±0.27 


4.83 ±0.17 


3.96 ± 0.29 


269 


0.10 ±0.01 


0.13 ±0.01 


1.09 ± 0.09 


5.49 ±0. 


,14 


4.43 ± 0.24 


1, 


.17 ±0.11 


3, 


.00 ±0.19 


7.16 ± 0.24 


4.85 ±0.19 


4.31 ±0.24 


292 






0.76 ±0.14 


4.55 ±0, 


,21 


3, 


.38 ± 0.36 





.70 ± 0.27 


2, 


.86 ±0.21 


4.52 ± 0.39 


4.67 ± 0.33 


2.02 ± 0.28 


298 


0.08 ± 0.01 


0.12 ±0.01 


1.24 ±0.12 


5.51 ±0, 


,16 


4, 


.33 ±0.29 


1, 


.21 ±0.12 


3.41 ±0.17 


6.97 ±0.81 


5.45 ± 0.20 


4.18 ±0.40 


300 


0.07 ±0.01 


0.12 ±0.01 


0.92 ±0.11 


5.30 ±0. 


,21 


4, 


.34 ± 0.35 


1, 


.43 ± 0.37 


3 


.15 ±0.19 


5.75 ± 0.36 


5.03 ±0.29 


4.20 ±0.27 


303 


0.04 ± 0.01 


0.06 ± 0.01 


0.92 ±0.20 


5.42 ±0, 


,20 


3 


.90 ± 0.62 


1, 


.07 ± 0.32 


3, 


.12 ±0.21 


5.15 ±0.70 


4.80 ±0.27 


3.32 ±0.37 


309 


0.08 ± 0.01 


0.12 ±0.01 


0.97 ±0.12 


5.61 ±0, 


,19 


3, 


.75 ±0.32 


0, 


,75 ± 0.25 


3, 


,27 ±0.18 


7.00 ± 0.32 


5.04 ± 0.24 


4.00 ± 0.27 


353 






0.75 ±0.12 


5.25 ±0, 


,16 


3, 


,95 ±0.26 


0, 


,88 ± 0.16 


3, 


,06 ± 0.16 


6.04 ±0.31 


4.43 ±0.21 


3.78 ± 0.20 


359 






0.85 ±0.14 


4.79 ±0, 


,21 


3, 


.73 ± 0.46 


1, 


,05 ± 0.32 


3, 


,25 ± 0.20 


5.36 ± 0.39 


3.84 ± 0.33 


3.40 ± 0.32 


360 






-0.15 ±0.66 


4.10 ±0, 


,86 


2, 


,95 ± 1.61 


1, 


,13 ± 0.92 


2, 


,40 ± 0.75 


2.44 ± 1.20 


8.18 ±0.89 


4.36 ±0.72 


375 


0.12 ±0.01 


0.15 ± 0.01 


0.95 ±0.10 


5.40 ± 0, 


,15 


4, 


,37 ±0.26 


1, 


,19 ± 0.09 


3, 


.48 ± 0.14 


6.82 ± 0.27 


5.50 ±0.41 


4.45 ± 0.37 


381 


0.04 ±0.03 


0.07 ±0.02 


0.91 ± 0.12 


5.15 ±0, 


,27 


4, 


,02 ±0.29 


0, 


.99 ± 0.18 


2, 


.98 ± 0.22 


5.00 ±0.50 


4.31 ±0.34 


3.87 ±0.29 


391 


0.05 ±0.01 


0.08 ± 0.01 


1.22 ±0.14 


5.07 ±0, 


,18 


4, 


,00 ±0.32 


1, 


,02 ± 0.24 


3, 


.15 ± 0.19 


5.41 ± 0.48 


4.50 ±0.36 


3.54 ±0.43 


408 


0.03 ±0.01 


0.06 ± 0.01 


0.46 ± 0.33 


5.66 ±0, 


,23 


3, 


.94 ± 0.32 


1, 


,10 ±0.28 


3, 


,19 ± 0.20 


6.20 ± 0.34 


4.18 ±0.33 


3.66 ±0.39 


409 


-0.01 ±0.01 


0.02 ± 0.02 


0.76 ±0.32 


4.88 ±0, 


,53 


2, 


,55 ±0.68 


0, 


,71 ± 0.40 


4, 


,19 ± 0.46 


4.07 ±0.92 


3.57 ±0.75 


4.24 ±0.70 


410 






1.08 ±0.21 


5.18 ±0, 


,.36 


3 


,62 ± 0.48 


1, 


,21 ± 0.29 


3, 


,93 ± 0.34 


6.75 ± 0.52 


4.35 ±0.41 


4.28 ±0.75 


412 


0.01 ±0.01 


0.03 ±0.01 


1.19 ±0.40 


5.32 ±0, 


,24 


3, 


.30 ± 0.32 


0, 


.95 ± 0.27 


2, 


.57 ± 0.23 


4.42 ± 0.52 


4.38 ±0.39 


3.35 ±0.24 


463 






1.02 ±0.11 


5.10 ±0. 


,18 


4, 


.16 ±0.28 


1, 


.09 ±0.17 


3, 


.31 ±0.19 


6.20 ± 0.75 


4.46 ± 0.24 


3.99 ±0.30 


481 


-0.01 ±0.01 


0.03 ± 0.01 


1.23 ±0.23 


4.49 ± 0. 


,42 


2, 


.53 ±0.71 


1, 


.14 ± 0.39 


2, 


.87 ± 0.42 


3.66 ± 0.75 


4.64 ± 0.60 


3.47 ± 0.45 


493 


-0.07 ±0.02 


-0.04 ± 0.02 


0.44 ± 0.34 


5.41 ±0. 


,56 


2, 


.28 ± 0.94 


1, 


.39 ± 0.55 


3, 


.06 ± 0.59 


0.88 ± 1.14 


3.99 ± 0.74 


3.59 ±1.07 


531 


0.08 ± 0.00 


0.12 ± 0.00 


0.99 ± 0.09 


5.33 ± 0. 


,14 


4, 


.54 ±0.18 


1, 


.16 ±0.09 


3, 


.15 ±0.14 


6.35 ± 0.22 


4.63 ± 0.20 


4.37 ± 0.23 


534 








4.15 ±0, 


,35 


3, 


.07 ±1.72 






2, 


.60 ± 0.38 


3.80 ± 0.92 


1.93 ±0.85 




536 


0.07 ±0.01 


0.10 ±0.01 


1.11 ±0.11 


5.30 ±0, 


,17 


3, 


.99 ±0.29 


1, 


.16 ±0.13 


3 


.06 ±0.16 


6.04 ± 0.77 


4.74 ± 0.23 


3.88 ± 0.28 
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TABLE A3 
First Partial Derivatives of the Models 



X 


ax 


ax 


ax 


ax 


ax 


ax 


a lost 


e[z/Hi 


S[a/Fo] 


a[a/N] 


a[a/c] 


S[a/Cal 


H<5f 


-1.813 


-1.356 


2.064 








H<5a 


-4.299 


-4.185 


4.743 








CNi 


0.122 


0.198 


0.058 


-0.147 






CNa 


0.126 


0.213 


0.066 


-0.164 


-0.507 




Ca4227 


1.000 


1.444 


-0.034 


0.421 


2.259 


-1.767 


G4300 


1.462 


1.748 


1.584 




-4.779 




H7F 


-2.612 


-2.232 


1.181 








H7A 


-4.229 


-4.637 


3.788 








Fe4383 


2.156 


4.920 


-4.488 








Ca4455 


0.693 


1.434 


0.104 








Fe4531 


0.973 


2.084 


-1.013 








C4668 


2.142 


8.953 


1.209 




-50.216 




H/3 


-1.154 


-0.629 


0.289 








Fe5015 


1.001 


3.208 


-1.574 








Mgb 


1.761 


3.183 


2.123 









Note. — The partial derivatives given above are computed by adding 0.3 dex to the adopted reference parameters 
{logt, [Z/H], [a/Fe], [a/N], [q/C], [a/Ca]} = {-0.15,0.3,0.2,0,0,0}. Note that we defi ne age in units of 10^° yr, and thus 
the reference age is equivalent to t = 7 Gyr, or a redshift of formation of zj = 2.4 (Kelso n et alJl2001^ . Thus dX/dp = 
AX/0.3, where p £ {logf, [Z/H], [a/Fe], fa/N], [a/C], [a/Ca]}. Where the partial derivatives are not given, and assumed 
to be identically zero, Thomas, Maraston, & Bonder (2003a) and Thomas. Maraston, & Korn (2004) had either reported, 
qualitatively, that the index was insensitive to changes in the relevant parameter, or the model cannot contain information on 
the sensitivity of the index to the parameter because they did not (or could not) include it in their analyses. This fact helps 
to underscore the importance of zero-pointing the mo del to a fiducial observed d ataset as it re moves the importance of such 
unkno wn /undefined dependences to zeroth order. The lThomas. Maraston. fc Ben der (2003^) and lThomas. Maraston. fc KornI 
ll2004fl set of predictions is highly non-linear such that at other locations in the 6D parameter domain these derivatives are 
not valid. As a result, the above partial derivatives should only illustrate a broad sensitivity of the indices to changes in the 
model parameters. 

TABLE A4 

The Zero-point Offsets for Recalibrating the Models to the Data 



X 


-^modcl 


-''^(CLiass) 


Ax 


units 




0.261 


0.519 


0.258 


A 


U5a 


-2.109 


-1.468 


0.641 


A 


CNi 


0.080 


0.060 


-0.020 


mag 


CNa 


0.113 


0.094 


-0.020 


mag 


Ca4227 


1.432 


0.958 


-0.474 


A 


G4300 


6.042 


5.326 


-0.716 


A 


H7F 


-1.461 


-1.265 


0.196 


A 


H7A 


-5.829 


-5.004 


0.825 


A 


Fe4383 


5.496 


4.138 


-1.358 


A 


Ca4455 


1.945 


1.096 


-0.848 


A 


Fe4531 


3.573 


3.166 


-0.407 


A 


C4668 


7.228 


6.074 


-1.154 


A 


H/3 


1.755 


1.729 


-0.025 


A 


Fe5015 


5.855 


4.690 


-1.165 


A 


Mgb 


4.422 


4.000 


-0.422 


A 



Note. — The model values given above are computed at the adopted reference parameters 

{logt, [Z/H], [a/Fe], [a/N], [a/C], [a/Ca]} = {-0.15,0.3,0.2,0,0,0}. Note that we define age in unit s of 10^" yr, and 
thus the reference age is equivalent to t = 7 Gyr, or a redshift of formation of Zf = 2.4 iKelson et alJl2001^ . 



TABLE A5 

Tested Variations in the Zero-point of the Models 



Model 


log[i/{10 Gyr)] 


[Z/H] 


[a/Fe] 


[a/N] 


[a/C] 


[a/Ca] 


Qa 


-0.15 


0.30 


0.20 


0.00 


0.00 


0.00 


A 


0.00 


0.30 


0.30 


0.00 


0.00 


0.00 


B 


-0.15 


0.30 


0.00 


0.00 


0.00 


0.00 


C 


-0.15 


0.00 


0.00 


0.00 


0.00 


0.00 



Note. — "The zero-point used in the analysis presented in the paper. 



